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A  Fortran  computer  program  to  calculate  radar  maximum  range,  written  for  the  NRL  CDC-3800 
computer  but  adaptable  to  any  computer  with  a  Fortran  compiler,  is  described.  The  computation  follows 
previously  established  principles,  with  the  pattern  propagation  factors  set  equal  to  one,  so  that  the  range 
calculated  Is  for  free  space  in  the  sense  that  earth’s  surface  effects  are  not  taken  Into  account.  However, 
the  effects  of  a  standard  atmosphere  are  included  in  the  calculation.  Reflection-Interference  effects  can  be 
separately  described  by  utiHeing  the  calculated  free-space  range  as  an  input  to  computer  plotting  programs. 

The  program  calculates  the  range  for  any  specified  probability  of  detection,  false-alarm  probability, 
and  Swerling  fluctuation  case  by  utilizing  a  slightly  modined  subroutine  written  by  Fehiner  and  coworkers 
of  the  Johns  Hopkins  University  Applied  Physics  Laboratory.  Postdetection  (noncoherent)  integration  is 
assumed.  The  system  noise  temperature  is  computed  including  effects  of  galactic,  cosmic  blackbody,  solar, 
and  tropospheric  noise,  and  the  tropospheric  rnolecular  absorption  for  oxygen  and  water  vapor  Is  calculated 
for  a  sUndard  atmosphere.  The  effect  of  refraction  on  the  ray  path  is  included  in  the  absorption  calculation 
by  ray  tracing,  assuming  a  negative-exponential  refractivity-height  profile.  The  range  of  validity  of  the  noise 
temperature  and  absorption  calculations  is  approximately  100  MHz  to  100  uHz.  The  computation  requires 
a  few  seconds  with  the  CDC-3800  computer. 
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ABSTRACT 


A  Fortran  computer  program  to  calculate  radar  maximum  range, 
written  for  the  NRL  CDC-3800  computer  but  adaptable  to  any  computer 
with  a  Fortran  compiler,  is  described.  The  coi:;puUition  follows  previously 
established  principles,  with  the  pattern-propagation  factors  set  equal  to  one , 
so  that  the  range  calculated  is  for  free  space  in  the  sense  that  earth’s  sur¬ 
face  effects  are  not  taken  into  account  However,  the  effects  of  a  standard 
atmosphere  are  included  in  the  calculation.  Reflection-interference  effects 
can  be  separately  described  by  utilizing  the  calculated  free-space  range  as 
an  input  to  computer  plotting  programs. 

The  program  calculates  the  range  for  any  specified  probability  of  detec¬ 
tion,  false-alarm  probability,  and  Swerling  fluctuation  case  by  utilizing  a 
slightly  modified  subroutine  written  by  Fehlner  and  coworkers  of  the 
Johns  Hopkins  University  Applied  Physics  Laboratory.  Postdetection 
(noncoherent)  integration  is  assumed.  The  system  noise  temperature  is 
computed  including  effects  of  galactic,  cosmic  blackbody,  solar,  and  tropo¬ 
spheric  noise,  and  the  tropospheric  molecular  absorption  for  oxygen  and 
water  vapor  is  calculated  for  a  standard  atmosphere.  The  effect  of  refrac¬ 
tion  on  the  ray  path  is  included  in  the  absorption  calculation  by  ray  trac¬ 
ing,  assuming  a  negative-exponential  refractivity-height  profile.  The  range 
of  validity  of  the  noise  temperature  and  absorption  calculations  is  approx¬ 
imately  100  MHz  to  100  GHz,  The  computation  requires  a  few  seconds 
with  the  CDC-3800  computer. 
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A  FORTRAN  COMPUTER  PROGRAM  TO  CALCULATE  THE  RANGE 

OF  A  PULSE  RADAR 


INTRODUCTION 

The  program  to  be  described  performs  a  calculation  of  the  free-space  range  of  a 
monostatic  pulse  radar,  following  the  principles  presented  in  NRL  Reports  6930  and 
7010  (1).  The  range  calculated  is  “free  space"  in  the  sense  that  the  effects  of  the  earth’s 
surface,  such  as  reflection-interference  and  below-the-horizon  shadow,  are  not  taken  into 
account.  However,  the  absorption  and  noise  of  a  standard  atmosphere,  galactic  noise, 
and  solar  noise  are  taken  into  account.  For  detailed  definitions  of  some  of  the  quantities 
involved  in  range  calculation  and  for  a  discussion  of  the  theory,  the  reader  is  referred  to 
those  reports.  In  NRL  Report  6930,  a  Range-Calculation  Work  Sheet  was  presented  to 
systematize  the  calculation  and  thus  to  simplify  handling  the  rather  large  number  of 
quantities  and  computational  steps  involved.  The  use  of  a  computer  program  represents 
a  still  further  step  toward  simplification  of  the  calculation,  and  also  minimizes  the  possibility 
of  error.  The  computation  requires  punching  one  card  with  the  input  data  (radar  parameters 
and  related  quantities).  To  guard  against  error  in  punching  this  card,  the  program  prints 
out  all  of  the  input  data  as  well  as  the  calculated  results.  The  program  has  been  given  the 
Fortran  name  RGCALC. 

Specifically,  this  program  is  a  Fortran  formulation  of  Eq.  (12)  of  NRL  Report  6930 
(Eq.  (3)  of  this  report),  with  the  pattern-propagation  factors  omitted.  The  signal  at  the 
target  is  therefore  assumed  to  be  due  .solely  to  direct-path  propagation,  and  if  the  specified 
transmitting  and  receiving  antenna  gains  are  those  of  the  beam  maxima,  the  target  is  assumed 
to  be  in  the  beam  maxima.  As  discussed  in  NRL  Report  6930,  the  equation  is  based  on 
the  assumption  that  the  detection  range  is  limited  by  the  normal  system  noise  —  i.e.,  that 
there  is  no  interference  from  manmade  signals  or  noise,  and  no  clutter  signals  caused  by 
echoes  from  extraneous  targets,  such  as  rough  sea  or  terrain,  rain,  or  any  profusion  of 
other  targets  in  the  vicinity  of  the  target  whos"  detection  is  being  considered.  The  range 
thus  calculated  may  be  called  the  “basic”  detection -range  capability  of  the  radar. 

The  Range-Calculation  Work  Sheet  of  NRL  Report  6930  requires  some  auxiliary 
calculations  and  the  use  of  some  sets  of  curves  to  determine  the  visibility  factor  (minimum- 
detectable  signal-to-noise  ratio);  the  antenna,  transmission-line,  receiver,  and  system  noise 
temperatures;  and  the  atmospheric  absorption  loss.  Calculations  equivalent  to  using  these 
curves  and  auxiliary  calculations  are  performed  within  FTogram  RGCALC.  The  only 
auxiliary  calculation  required  is  that  of  number  of  pulses  on  target,  for  a  scanning  radar. 

It  was  not  considered  feasible  to  do  this  calculation  in  the  computer  program  because  the 
number  of  pulses  is  sometimes  determined  by  a  signal  processor  rather  than  by  the  scanning. 
Hovrever,  the  calculation  is  not  difficult,  in  the  scanning  radar  case,  as  will  be  discussed 
later  in  this  report.  Postdeteciion  (noncoherent)  integration  of  the  pulses  is  assumed.* 


•See  Ref.  1  (NRL  Report  6930),  p.  18. 
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The  computer  program  can  also  be  adapted  to  calculating  the  range  of  a  CW  radar 
and  of  a  bistatic  radar  by  suitable  redefinition  of  some  of  the  input  parameters.  For  CW 
radar,  the  pulse  length  t  can  be  interpreted  as  the  effective  sampling  time  of  the  detection 
processor  m  nucroseconds;  or  if  such  a  processor  is  not  used,  the  parameter  ?  can  be  in- 
lerpreied  as  Lht  ivciprocal  of  tiiC  it-ceiver  prcdctcction  bandwidth  in  meiJahertz,  or  as  the 
length  of  time  a  scanning  beam  remains  on  the  target,  whichever  time  is  the  shorter.  The 
bandwidth  correction  factor  can  be  used  as  a  correction  for  non-optimum  processing  or 
filtering  of  the  signal.  The  “number  of  pulses  integrated”  should  be  set  equal  to  1.  The 
transmitter  power  Pf  is  defined  for  this  calculation  as  the  average  transmitted  power  in 
kilowatts. 

If  the  radar  is  bistatic,  the  only  reinterpretation  required  is  in  the  meaning  ui  the 
calculated  range  number.  Instead  of  representing  the  monostatic  range,  it  represents  the 
square  root  of  the  pnxluct  of  the  transmitter-to-target  range  Rf  and  the  target -to-recoiver 
range  R^.  In  other  words,  the  calculated  number  is  the  geometric  mean  of  i?,  and 
The  transmit  and  receive  antenna  gains  are  already  specified  separately  in  the  program 
because  they  are  actually  different  even  for  some  monostatic  radars.  The  target  cross 
section  a  must  of  course  be  the  bistatic  value. 

The  frequency  range  over  which  the  program  may  be  used  is  from  about  100  MHz 
to  100  GHz,  but  it  could  be  extended  downward  to  about  30  MHz  and  upward  to  perhaps 
150  GHz  without  incurring  gross  errors.  (Below  30  MHz  the  occurrence  of  ionospheric 
effects  and  above  150  GHz  the  multiplicity  of  water-vapor  absorption  resonances  invalidate 
the  equations  used  in  the  program.) 

The  program  is  written  with  the  option  of  calculating  range  for  either  specified 
detection  and  false-alarm  probabilities,  or  for  a  specified  signal-to-noise  power  ratio  (expressed 
in  decibels).  The  latter  option  is  useful  in  calculating  the  maximum  range  of  a  tracking 
radar  (as  distinct  from  a  search  or  acquisition  radar)  when  the  minimum  signal-to-noise 
ratio  for  successful  tracking  is  known. 

Except  for  the  fact  that  the  external  noise  from  celestial  and  terrestrial  sources  and 
the  absorpiion  that  occurs  in  the  earth’s  atmosphere  are  taken  into  account,  the  free  space 
range  of  the  radar  is  calculated.  Actually  it  would  be  more  accurate  to  call  this  range  a 
quasi-free-space  range  because  of  the  inclusion  of  celestial  and  terrestrial  noise  and  absorp¬ 
tion  effects.  Because  the  absorption  by  the  troposphere  is  dependent  on  the  elevation 
angle  of  the  ray  path,  the  target  elevation  angle  is  one  of  the  input  quantities  for  the 
range  calculation.  The  radar  is  assumed  to  be  located  at  or  near  the  earth’s  surface  — 
within  say  a  thousand  feet  of  sea  level.  Range  calculations  applicable  to  extraterrestrial 
locations  (e.g.,  satellite  or  space-ship  radars)  tan  be  made  by  setting  the  elevation  angle  to 
a  high  value,  e.g.,  90  degrees.  At  this  elevation  angle,  the  absorption  is  usually  negligible 
for  frequencies  appreciably  below  the  22-GHz  water-vapor  resonance  line.  Also,  the 
computed  absorption,  in  decibels,  is  printed  out,  so  that  the  computed  range  can  be 
revised  to  correct  for  it  if  a  true  free-space  range  is  desired.  The  ray  path  for  the  absorption 
calculation  is  computed  by  a  ray-tracing  algorithm,  assuming  an  exponential  decrease  of 
the  refractivity  with  height  (CRPL  Exponential  Atmosphere),  with  surface  lefractivity  of 
313  N  units  (2). 

The  principal  non-free-space  factor  not  taken  into  account  is  the  effect  of  the  reflec- 
t-on  and  absorption  by  the  earth’s  surface.  These  effects  may  modify  the  free-space  range 
greatly  under  some  conditions,  but  they  cannot  be  readily  taken  into  account  by  a  single 
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range  calculation.  To  depict  the  non-free  space  performance  of  a  radar,  what  is  needed  is 
a  graphical  representation  of  the  detection -range  contours  in  a  vertical  plane,  or,  for  a 
target  at  constant  altitude,  a  plot  of  the  variation  of  signal  strength  with  range,  relative 
to  the  signal  required  for  detection. 

Program.s  to  produce  such  plots  have  been  written,  and  have  been  described  in  a 
separate  report  (3).  One  of  the  inputs  required  for  these  programs  is  the  free-space  range 
of  the  radar.  Hence,  the  present  program  supplements  the  plotting  programs.  They  could 
be  combined  into  a  single  package.  However,  this  has  n>'i,  been  done  because  ordinarily 
it  is  not  objectionable  to  calculate  the  free-space  performance  and  to  plot  the  non-free- 
space  detection  curves  as  two  successive  operations. 

For  radars  whose  antenna  pattern  is  a  narrow  elevatable  beam,  no  significant  reflection- 
interference  effects  occur  when  the  beam  is  elevated  by  one  beamwidth  or  more,  and  the 
quasi-free-space  range  calculation  then  applies  directly. 

The  program  to  be  described  has  been  used  and  extensively  tested  over  a  period  of 
time.  Results  agree  with  manual  calculations  using  the  range-calculation  worksheet  of 
NRL  Report  6930  (1).  The  program  was  used  to  compute  the  ranges  of  actual  Navy 
radars  for  a  forthcoming  NRL  Radar  Division  report."'  The  execution  time  for  a  single 
radar  range  calculation,  for  all  five  Swerling  fluctuation  cases,  is  approximately  2  sec  on  the 
NRL  CDC-3800  computer,t  not  including  compilation  time  which  is  about  100  sec. 
(Compilation  can  of  course  be  avoided  by  having  an  “object  deck”  or  machine-language 
deck  punched,  and  using  it  instead  of  the  Fortran  or  “soiuce  deck.”) 


COMPARISON  WITH  OTHER  PROGRAMS 

Some  years  ago,  a  computer  program  was  developed  for  the  calculation  of  maximum 
radar  rays  by  a  contractor!  for  the  Scientific  and  Technical  Intelligence  Center  of  the 
Office  of  Naval  Intelligence  (ONI-STIC-50).  The  work  was  completed  about  1966.  This 
program  was  based  on  NRL  Report  5868,  an  earlier  edition  of  NRL  Report  6930  (1). 

It  utilized  some  ciu^es  published  in  that  report  by  reading  a  finite  number  of  data  points 
into  the  computer  and  interpolating  between  them;  this  was  done  for  the  “visibility 
factor”  (minimum-detectable  signal-to-noise  ratio),  the  antenna  noise  temperature,  and 
the  atmospheric  absorption  loss.  Decause  of  this  the  program  was  limited  to  calculating 
the  range  for  0.5  probability  of  detection  for  a  nonfluctuating  target,  and  to  the  frequency 
range  100  MHz  to  10  GHz. 

Another  computer  program  has  been  described  by  Boothe  (4).  This  program 
computes  the  probability  of  detection  as  a  function  of  the  range,  rather  than  computing 


•This  will  be  the  fifth  editiun  of  NRL  Report  5CJT,  4lh  ed.,  June  21.  1961.  "Navi'  Ra.lar  Systems 
Sur'-ey,"  R.D.  Tompkins. 

tThe  execution  time  depends  partly  on  the  number  of  pulses  integrated  end  on  whether  oi  not  other 
immediately  preceding  calculations  have  been  made  for  the  same  number  of  pulses  and  the  same  fals.  - 
alarm  probability. 

tTlie  contractor  was  Control  Data  Corporation,  and  the  programming  was  done  by  Mrs.  Irma  Wacbtel. 
Consulting  assistance  on  the  radar-equation  aspects  of  the  work  was  furnished  by  the  author  of  this 
report.  The  project  was  initiated  arvd  monitored  by  LCDR  William  Barron  of  ONI-STIC.  The  program 
IS  described  in  an  internal  ONI  document  5TiC-CW' O.o-l -66  titied  "Compuicr  Program  R-nO,  Radar 
Maximum  Range  Calculation,"  Feb.  23,  1967. 
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the  range  for  a  specified  probability.  Howcvci,  the  nraximum  range  for  specified  probability 
is  found  by  computing  probability  for  decreasing  range  values  until  the  desired  probability 
is  reached.  The  program  also  makes  use  of  atmospheric  absorption  curves  from  NRL 
Report  5868,  entered  as  data  into  the  computer;  consequently  it  is  limited  to  the  frequency 
range  100  MHz  to  10  GHz,  plus  perhaps  a  few  “spot’'  frequencies  up  to  100  GHz.  The 
report  does  not  state  the  method  used  for  evaluating  antenna  noise  temperature,  and  a 
program  listing  is  not  given.  The  effect  of  target  aspect  variation  on  cross  section  is  taken 
into  account  deterministically,  rather  than  statistically  using  Swerling’s  theory.  The  signal- 
to-noise  ratio  and  the  resulting  probability  of  detection  are  calculated  at  ranges  that 
decrease  in  steps  corresponding  to  observation  of  a  target  approaching  the  radar.  The 
target  is  assumed  to  be  changing  aspect  according  to  some  known  prescription,  during  this 
approach,  and  the  corresponding  cross-section  variation  is  calculated.  (A  missile  target  is 
assumed  in  Boothe’s  analysis.)  As  the  target  approaches,  when  the  probability  reaches  the 
specified  value,  the  range  is  printed  out  or  otherwise  recorded.  Either  single-scan  or 
cumulative  probability  can  be  specif ieil. 

D.  M.  White  has  described  a  comprehensive  program  (5)  to  analyze  radar  performance 
in  a  dynamic  situation,  computing  signal-to-noise  ratio  and  detection  probability  as  a  func¬ 
tion  of  time  and  target  position,  taking  into  account  the  effects  of  multipath  interference, 
clutter  echoes  from  the  sea  or  ram,  and  target  motion.  In  short,  the  program  simulates 
in  as  much  detail  as  is  practical  the  complete  radar-target  engagement,  for  a  single  target. 
This  program  utilizes  a  subroutine  written  by  L.  F.  Kehlner'*'  of  the  Johns  Hopkins  Applied 
Physics  Laboratory,  and  described  by  Fehlner  in  a  previous  report  (6),  to  calculate  the 
probability  of  detection  for  a  specified  signal-to-noisc  ratio,  false-alarm  number,  number 
of  pulses  integrated,  and  target  fluctuation  characteristic.  Any  one  of  five  fluctuation 
cases  can  be  specified;  the  nonfluctuating  case  and  the  four  Swerling  fluctuation  cases  (1). 
White  mentions  other  programs  that  have  been  written  by  Kirkwood  (7)  and  by  Nolen  (8). 

Killinger  (9)  has  developed  a  computer  program  that  calculates  the  ratio  of  signal  to 
noise-plus-clutter  as  a  function  of  t.arget  r-nge.  Probabilities  of  detection  and  false  alarm 
are  also  computed.  Maximum  detection  range  can  be  found  for  a  specified  signal  to  noise- 
plus-clutter  ratio. 

The  philosophy  of  the  program  to  be  described  in  this  report  is  somewhat  different 
from  those  discussed  above.  It  is  not  intended  to  simulate  the  radeir  performance  in  a 
dynamic  situation.  Instead,  it  is  intended  to  provide,  for  a  specified  target  size,  fluctuation 
model,  and  detection  probability,  a  single  number  that  will  serve  as  an  index  of  the  radar’s 
range  capability  —  a  “figure  of  merit.”  The  geophysical  enenronment  is  taken  into  account 
as  realistically  as  possible  except  that  effects  of  clutter,  rain,  and  multipath  interference 
are  omitted.  The  factors  that  are  believed  to  be  more  realistically  or  accurately  calculated 
than  in  other  programs  are  the  system  noise  temperature  (or  more  specifically,  the  tropo¬ 
spheric,  solar,  and  galactic  contributions  to  the  antenna  temperature)  eind  the  tropospheric- 
absorption  loss  (due  to  collision-broadened  absorption  resonances  of  the  oxygen  and  water- 
vapor  molecules).  The  antenna  noise  temperature  and  the  atmospheric  attenuation  arc 
computed  directly  rather  than  by  interpolation  using  data  entered  from  precalculated 
curves  or  tables;  consequently,  the  permissible  range  of  frequency  is  much  greater  than 
for  most  programs  using  precalculated  temperature  ^nd  absorption  data. 


Assisted  by  R.G.  Roll  and  G.T.  'rVuttor. 
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The  program  to  be  described  in  this  report  computes  the  detection  range  for  either 
a  steady  (nonfluctuating)  target  or  for  any  of  the  four  Swerling  fluctuation  cases  (or  for 
all  five  cases),  for  a  specified  probability  of  detection  and  a  specified  false-alarm  probability. 
Fehlner’s  subroutine  (which  he  named  MARCUM)  has  been  incorporated  into  the  NRL 
program,  with  slight  modifications,  for  this  purpose.  The  principal  modification  has  been 
to  provide  for  calculating  on  the  basis  of  false-alarm  probability,  rather  than  Marcum’s 
false-alarm  number.  Another  modification  insures  that  when  successive  calls  are  made  to 
the  subroutine  with  the  same  false-alarm  and  number -of-pulses  parameters,  the  bias-level 
calculation  is  not  repeated.  This  saves  an  appreciable  amount  of  computing  time  in  the 
iterative  procedure  used  to  determine  signat-to-noise  ratio  for  specified  probability.  (Sub¬ 
routine  MARCUM  actually  does  the  inverse  problem  of  computing  probability  for  a  specified 
signal-to-noise  n^tio.)  Because  of  these  and  other  changes,  the  subroutine  as  actually  used 
in  the  NRL  program  has  been  renamed  MARSWR  (acronym  for  Marcum-Swerling);  but  it 
is  basically  Fehlner’s  MARCUM  subroutine.  The  calculation  is  made  assuming  a  square- 
law  detector,  whereas  most  radar  receivers  employ  a  linear-rectifier  detector,*  but  the 
difference  in  performance  oi  the  two  detector  types  is  about  0.2  dB  at  most,  depending 
on  number  of  pulses  integrated. 


DATA  INP'JTS 

The  inputs  to  Program  RGCALC  for  a  single  radar  range  calculation  are  punched  on 
a  single  data  card.  The  format  specifications  for  this  card  are  as  follows.  Each  of  the 
listed  quantities  is  discussed  briefly  in  the  following  paragraphs. 


Data  Item 

Format  Specification 

Card  Columns 

Transmitter  power,  kW 

F6.0 

1-6 

Pulse  length,  ^sec 

F6.0 

7-12 

Transmit  antenna  gain,  dB 

F4.0 

13-16 

Receive  antenna  gain,  dB 

F4.0 

17-20 

Target  cross  section, 

F6.0 

21-26 

Frequency,  MHz 

F6.0 

27-32 

-Antenna  ohmic  loss,  dB 

F4.0 

33-36 

Receiving  line  loss,  dB 

F4.0 

37-40 

Transmit  line  loss,  dB 

F4.0 

41-44 

Antenna  pattern  scan  loss,  dB 

F4.0 

45-48 

Miscellaneous  loss,  dB 

F4.0 

49-52 

Bandwidth  factor,  dB 

F4.0 

53-56 

Receiver  noise  factor,  dB 

F4.0 

57-60 

Number  of  pulses 

15 

61-65 

Probability  of  detection 

F4.J 

66-69 

False-alarm  exponent 

F4.0 

70-73 

Swerling  fluctuation  case 

11 

14 

Target  elevation  angle,  deg 

F4.0 

7 ',-7  8 

Galactic  noise  code 

12 

79-80 

•See  Ref.  1  (NRL  Report  69.70),  p.  29. 
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The  data  items  of  the  preceding  list  denoted  “probability  of  detection"  (Cols.  G6-69) 
and  “false-alarm  exponent”  (Cols.  70  73)  are  actually  so  defined  only  if  the  “case"  parameter 
(Col.  74)  is  0,  1,  2,  3,  4,  or  5.  If  the  rase  parameter  is  6  c.  7,  as  will  be  discussed  sub¬ 
sequently  in  more  detail,  the  range  is  calculated  for  a  specified  signal-to-noise  ratio  in 
decibels  (Case  6)  or  ratios  (Case  7).  If  one  signal-to  noise  ratio  is  to  be  specified,  it  goes 
in  card  columns  66-69  in  place  of  the  probability  of  detection.  The  second  signal-to-noise 
ratio  goes  in  Cols.  70-73,  in  place  of  the  false-alarm  exponent. 

Data  items  shown  as  having  an  F  format  specification  can  be  entered  as  a  number 
including  a  decimal  point,  with  the  number  positioned  anywhere  within  the  card-column 
field.  If  an  F-specificatiou  number  happens  to  be  an  integer,  it  can  also  be  entered  without 
a  decimal  point,  but  in  that  case  i:  must  be  right-adjusted  within  the  card-column  field. 

The  specification  f6.0  means  a  number  of  6  characters  or  less  including  decimal  jxrint  and 
sign,  if  any  (positive  sign  is  implied  if  no  sign  is  given).  The  decimal  point  can  be  positioned 
anywhere  in  the  field.  Data  items  having  an  I  format  specification  must  be  integers  (no 
decimal  point),  and  must  be  right-adjusted  in  the  colrmn  field. 

This  single  card  contains  all  the  numerical  data  required  for  a  radar  range  calculation. 
However,  two  data  cards  are  required  for  each  radar  calculation;  the  other  (first)  card 
contains  any  alphanumeric  material  that  may  be  required  to  identify  the  radar.  This 
material  is  punched  anywhere  in  the  80-column  fie'd  of  the  data  card,  and  it  will  be  printed 
out  at  the  top  of  the  page  preceding  the  listing  of  input-output  quantities. 

Calculations  for  any  desired  number  of  radars  can  be  made  in  one  computer  run  by 
stacking  the  data  cards  in  the  following  manner: 

Cards  1,  3,  5,  7,  .  .  .  Alphanumeric  material  identifying  the  radars 

Cards  2,  4,  6,  8,  .  .  .  Data  cards  giving  numerical  parameters  corresponding  to  the 
preceding  identifier  cards. 

If  it  is  not  desired  to  provide  alphanumeric  identifying  material,  blank  cards  should  be 
inserted  at  positions  1,  3,  5,  7,  .  .  .  of  the  data  deck.  The  last  card  of  the  deck  is  an 
end-of-file  card.  When  this  is  encountered,  the  job  will  be  terminated. 

In  the  discussion  of  definitions  that  follows,  a  basic  principle  should  be  kept  in  mind. 

In  any  radar  system,  the  partitioning  of  the  system  into  sections  called  “antenna,”  “trans¬ 
mission  line,”  "receiver,”  and  “transmitter”  is  somewhat  arbitrary  (sec  NRL  Report  6930, 
p.  47,  Fig.  10,  and  NRL  Report  7010,  p  38,  Fig.  5).  The  points  in  the  system  at  which 
I '.is  arbitrary  partitioning  is  done  determines  the  numerical  values  of  losses,  gains,  power, 
and  noise  temperatures  to  be  assigned  to  the  factors  which  will  subsequently  be  identified 
as  La,  Lr,  L(,  P(,  G/,  G,,  Tq,  Tr,  and  Tg  {NF).  The  range  calculation  will  come  out  the 
same  no  matter  how  this  partitioi'ing  is  done  if  the  assignment  of  values  to  all  these 
quantities  is  consistent  with  the  partitioning  selected.  Values  of  losses  in  decibels  are  to 
be  entered  on  the  data  card  as  positive  numbers. 

Transmitter  Power.  See  NRL  Report  6930  (1),  p.  11.  Symbol  P/.  This  is  the  pulse 
power  of  the  radar  in  kilowatts. 


Piilso  Length.  So<'  NRl.  Rcporl  fiR-iO,  jv  11.  Tho  .•ivnilml  r  is  usi'd  for  tiu'  duraiion 
liolwpon  3-db  points  of  tIu'  trnnsniittod  RK  pul.e<\  in  mitrosci'oiui.s.  If  tIu'  radar  is  of  tho 
inilsoTOtnprossion  tyjK'.  tin’  urnompressfd  puls<-  It-nplli  applies,  assuniinu  that  P/  is  tin- 
power  of  tlio  uneonipressod  puls''.  Tho  hasie  rule  is  tnat  the  prodiu  t  I'/t  must  equal  the 
tran.smittwl  pulse  energy.  (Mori  spoeifieally,  the  pul.se  )  iwer  in  kilowatts  tinu’s  10-^ 
multiplied  by  the  pulse  length  in  mierosreonds  times  !0~6,  must  equal  the  transmitted 
pulse  energy  in  watt -.second.) 

.^ntenna  Gain.  Sio  N'RL  Report  6930  p.  12.  The  power  g.nns  of  the  transmitting 
antenna  [C/)  and  receiving  antenna  (G^)  are  in  deciliels.  Power  gain  G  is  to  lu’  distinguished 
from  directive  gain  I),  ihese  quantities  are  related  by  (.‘  -  !;D,  where  h  is  the  radiation 
efficiency  (/;  <1  1).  The  radiation  efficiency  is  a  measure  of  ohmic  or  heat  loss  in  the 
antenna,  and  should  not  be  confu.sed  with  aperture  efficiency,  which  measures  the  relatii'ii- 
ship  between  the  directive  gain  actually  obtained  and  that  which  would  have  1  e.’ii  obtained 
if  the  aperture  were  uniformly  illumii  ilod. 

Target  Cross  Section,  bee  NRL  Rei'ort  6930,  p.  13.  The  symbol  •'  is  useii  for  tho 
radaj  cross  section  of  the  targi’t  in  square  meters.  For  eomi'iu-ison  of  the  porformanee 
of  competing  sy.'tems,  the  value  <i  =  1  ni2  i.s  often  used. 

Frequency  See  NRL  Report  6930,  p.  14.  The  symbol  is  for  the  radar  frequency 
in  megdhetU. 

Antenna  Ohmic  Loss.  See  .N'RL  Report  6930,  p.  48.  The  symbol  is  for  the  ohmie 
loss  of  the  antenna  ex()ressed  in  decibels.  Kven  though  tliis  loss  is  taken  into  account  hy 
the  fact  that  G  represc  its  the  power  gain  rather  than  the  directive  gain  of  the  antenna,  it 
must  also  be  entered  separately  because  of  ii.s  contributioti  to  thi'  system  noise.  (Its 
inclusion  in  the  power  gain  accounts  only  for  its  effect  on  the  transmitted  and  received 
signal  powers.)  If  there  are  separate  transmitting  and  receivang  antennas.  Lq  refers  to  the 
receiving  antenna  only.  This  quantity  is  negligible  for  many  types  of  aiuennas,  particularly 
for  parabolic  reflector  types,  for  which  tne  appro .\imat ion  =  0  dB  is  usually  justifiable. 
Certain  types  of  airay  antennas,  especially  those  that  employ  fiequeiiey  or  jihase  scanning, 
may  have  appreciable  ohrreie  loss. 

Receiving  Line  Loss.  See  NRL  Report  6930,  p.  70.  The  syrr.liol  hr  is  used  for  the 
loss  of  the  receiving  transmission  line  id  decibels.  This  loss  usually  includes  duplexor  or 
circulator  losses;  the  prefatory  remark  concerning  partitioning  of  the  receiving  system 
applies. 

Transmitting  Line  Loss.  See  NRL  Rcporl  6930,  p.  70.  The  symbol  is  for  llie 
loss  of  the  transmitting  portion  of  the  trunsmi.ssion  line  in  decibels  (not  usually  identical 
to  Lr).  Duplexer  loss  is  usually  included.  I'he  remark  concerning  partitioning  of  the 
system  applies. 

Antenna-Patte'';i  Scan  Loss.  See  NRL  Report  6930,  p.  70.  Symbfd  This  loss 
reflects  the  facts  that  (a)  the  number  of  pulses  integrated  for  a  scanning  radar  is  somewlmt 
arbitrarily  taken  to  be  tlie  num'ner  occurring  while  the  target  is  within  the  half-i>ower 
beamwidth  of  the  antenna,  and  (b)  th?  beam  does  not  have  full  uniform  gain  within  this 
beamwidth  and  zero  gain  elsewhere.  For  a  non.scanning  radar  aimed  directlv  at  a  target, 
hp  =  0  dB.  For  a  simple  azimuth-scanning  radar,  hp  1 .6  dB.  For  a  simultaneously 
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azimuth-  and  pk'valion-scanmn^'  radai,  Lp  '  3.2  dB  i.s  a  reasonaiile  assumption,  although 
this  result  is  based  on  a  crude  rather  than  a  sophisticated  analysis. 

Miscellaneous  Loss.  See  NRL  Report  6930,  p.  82-84.  Symbol  L.v.  Among  the 
[lossible  losses  that  may  be  included  here  are  collapsing  loss,  signal-processing  loss,  array- 
fill-time  loss,  Ix-am-squint  loss,  [xilanzation-rotation  loss,  and  rain-absorption  loss  (if  the 
rainstorm  extent  is  less  than  the  radar-to-taiget  range).  The  decibel  value  of  this  loss  is 
obtained  by  directly  adding  the  decibel  values  of  individual  contributing  losses. 

Bandwidth  Factor.  See  WRL  Report  6930,  p.  14.  Symbol  Cg  is  for  the  decibel  loss 
resulting  from  a  mismatch,  in  the  North-filter  sense,  between  the  fiulse  characteristics  and 
the  receiver  filter  transfer  characteristic.  For  a  simple  pulse  radar,  this  relationship  can 
be  analyzed  adequately  in  terms  of  the  pulse  length  and  shape,  and  the  biter  bandwidth. 

For  most  radars  of  this  type  it  is  reasonable  to  assume  =  0  dB,  in  the  absence  of 
specific  knowledge  to  the  contrary.  For  pulse-compression  radars,  there  is  usually  some 
loss  associated  with  the  compression  filter,  ranging  from  perhaps  0.5  dB  to  several  decibels, 
defiending  on  the  technique  employed  and  the  compression  ratio.  As  is  done  with  the 
loss  factors,  Cf^  is  to  be  entered  as  a  positive  decibel  number. 

Receiver  Noise  Factor.  (.Also  called  receiver  noise  figure;  although  “fiev re”  is  perhaps 
the  most  common  usage,  lEEK  stanaards  give  preference  to  "factor.”)  See  NRL  Report  6930, 
p.  50.  Symbol  .VF  or  F„.  The  receiver  noise  factor  and  receiver  noise  temperature  are 
alternative  ways  of  e.xpressing  the  same  property  of  the  receiver,  but  the  noise  factor  has 
been  chosen  here  because  it  is  the  quantity  more  commonly  given  in  receiver  specifications. 
The  decibel  value  of  the  noi.se  factor  is  to  be  entered  on  the  data  card. 


Number  of  Pulses.  See  NRL  Report  6930,  pp.  71  and  72.  Symbol  A/.  If  this  number 
IS  determined  by  a  signal  processor,  it  must  be  found  by  reference  to  the  characteristics 
of  the  processor.  When  it  is  determined  by  the  scanning  action  of  the  radar  antenna,  and 
if  a  simple  azimuth  scan  is  employed,  the  appropriate  formula  is 

6  RPM  •  cos  b/ 


where  o  is  the  azimuthal  half-power  beamwidth,  degrees:  PRF  is  the  pulse  repetition 
frecjuency,  hertz;  RPM  is  the  rotation  rate  of  the  antenna,  revolutions  per  minute;  and 
d,.  IS  the  elevation  angle  of  the  target.  (The  term  cos  Op  is  significant  only  when  a  target 
IS  ,  ?n  elevation  angle  of  about  10  degrees  or  more.  For  vertical-fan -beam  radars  the 
.  g*’  IS  usually  calculated  at  an  elevation  angle  low  enough  so  that  cos  =  1.) 


For  simultaneous  azimuth-  and  elevation-scanning  radars,  assuming  a  sawtooth-motion 
elevation  scan  and  a  uniform-speed-rotation  azimuth  scan,  the  appropriate  formula  is 


<■■0  PRF 


M  =  ;r - 


6  u>p  •  Ip  ■  RPM  ■  cos  Op' 


(2) 


)!'  which  o,  PRF,  RPM,  and  Hp  have  the  same  definitions  as  before,  0  is  the  vertical  beam- 
width,  oj,.  is  the  vertical  .semning  speed  in  degrees  per  second  at  the  target  elevation  angle, 
and  t,  IS  the  veriicai-scan  period  in  seconds,  including  the  dead  time  it  any. 
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The  number  of  pulses  to  be  used  for  radars  of  other  scan  types  must  t)e  analyzed  on 
an  individual  basis,  as  discussed  in  NRL  Report  6930. 

The  number  of  pulses  is  entered  on  the  data  card  as  an  integer  (I-format).  Consequently 
the  number  must  be  right-adjusted  in  the  field  of  columns  61  through  65.  If  the  number 
calculated  from  the  above  formulas  is  not  an  integer,  it  should  be  rounded  off  to  the 
nearest  integer.  (The  subroutine  that  calculates  detection  probability  requires  an  integer 
for  the  number  of  pulses  integrated.) 

Probability  of  Detection.  See  NRL  Report  6930,  pp.  18  and  19.  Symbol  P^-  Prob¬ 
ability  is  here  given  in  the  mathematical  sense  of  a  number  between  0  and  1  (not  as  a 
percentage  figure).  Values  larger  than  0.99  should  not  be  entered  because  computational 
difficulties  result.  Also,  values  smaller  than  the  false-alarm  probability  are  meaningless; 
for  practical  purposes,  P^  should  be  at  least  an  order  of  magnitude  larger  than  Pfa-  (Ordi¬ 
narily  it  is  many  orders  of  magnitude  larger.)  Tyjiical  values  of  of  interest  range  from 
about  0.1  to  0.95. 

False-Alarm  Exponent.  See  NRL  Report  6930,  pp.  18-19.  Symbol  —  logjQ  Pfg. 

Typical  values  of  false-alarm  probability  range  from  10'  to  10'*^.  The  number  to  be 
entered  on  the  data  card  is  the  positive  value  of  the  exponent  (power  of  ten).  Thus,  for 
Pfg  =  10~®,  enter  the  number  6.0  on  the  data  card;  for  Pf^  =  2.5  X  10“®,  enter  5.6. 

Swerling  Fluctuation  Case.  See  NRL  Re))ort  6930,  p.  28.  Five  cases  are  considered, 
with  0  representing  the  nonfluctuating  target  and  integers  1  through  4  representing  the 
4  Swerling  cases  as  defined  in  NRL  Rejjort  6930  (and  elsewhere).  If  the  numbers  0,  1.2, 

3,  or  4  are  punched  on  the  data  card  in  Column  74,  the  corresponding  fluctuation  case 
will  be  calculated.  If  5  is  punched,  all  5  cases  will  be  calculated. 

Further  options  are  provided  by  a  6  or  a  7  punch  in  Column  74.  A  6  punch  signifies 
that  the  range  is  to  be  calculated  for  a  specified  signal-to-noise  ratio  rather  than  for 
specified  probabilities  of  detection  and  false  alarm.  For  this  case,  the  signal-to-noisc  ratio, 
in  decibels,  is  punched  in  Column  66-69,  where  probability  of  detection  would  ordinarily 
appear.  If  a  7  is  punched  in  Column  74,  the  calculation  of  range  will  be  made  for  two 
different  signal-to-noise  ratios,  one  given  in  Columns  66-69,  the  other  in  Columns  70-73. 

When  either  6  or  7  is  punched  in  Column  74,  the  number-of -pulses  entry.  Columns  61-65, 
is  ignored.  Likewise,  if  6  is  jiunched,  the  false-alarm  exponent  entry,  Columns  70-73,  is 
Ignored. 

Target  .,  /ation  Angle.  See  NRL  Report  6930,  pp.  48,  68,  69,  and  72  through  80. 

As  mentioned  in  the  Introduction,  this  factor  enters  into  the  “quasi-free-space”  range  cal¬ 
culation  because  the  effect  of  the  earth's  atmosphere  on  the  antenna  noise  temperature 
and  on  absorption  loss  is  taken  into  account.  The  elevation  angle  is  to  be  entered  in  degrees. 
If  a  range  calculation  applicable  in  empty  space  is  desired,  a  close  approximation  can  l>e 
obtained,  except  at  frequencies  near  the  water  vapor  and. near  or  above  the  oxygen  resonances 
(22  and  60  GHz,  respectively)  by  setting  the  elevation  angle  to  90  degrees,  because  for 
this  setting  the  ab.sorption  is  virtually  negligible.  Also,  since  the  calculated  absorption  is 
one  of  the  printed-out  quar  tities.  correction  for  it  can  )>e  made. 

Galactic  Noise  Code.  Kcc  NRL  Rei«jii  6930,  |j.  49,  Fig.  11.  As  shown  in  the  referenced 
figure,  the  noise  received  from  the  galaxy  to  w  hich  the  solar  system  belongs  varies  rlepending 
on  the  fiart  of  the  galaxy  toward  which  the  antenna  is  pointed.  This  direction  is  not 
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usually  predictahle.  Therefore  the  options  of  caleulatinK  the  ratlar  rarine  for  three  choiees 
of  galaetic  anJ  solar  noise  levels  are  provided.  The  <‘odt  s  are  -]  for  minimum  galactic- 
noise,  0  for  average  noise,  and  -t-l  for  maximum  '■.alaetic  nui.so.  (The  maximum  and  min- 
values  are  shown  by  dashed  line's  in  the-  referenc'txl  figure  of  NRL  Report  6930.)  The 
number  entered  is  to  be  right -adjusted  in  Columns  79-80 


PROGRAM  OUTPLT 

The  output  of  Program  RGCALC  is  a  single  printed  page  for  each  set  of  data  inputs 
(two  data  cards).  The  alphanumeric  material  of  the  first  data  card  is  printed  at  the  top 
of  the  page.  Then  the  numerical  input  data  are  printed,  both  as  a  record  of  the  data  and 
to  ensure  that  the  data  card  was  correctly  [lunched.  Next  arc'  printed  some  intermediate 
output  quantities  such  as  the  computed  noise  tem[ieratures  of  the  system  conqxments  and 
of  the  overall  receiving  system,  and  the  tro|X)spheric  absorption  for  a  two-way  path  through 
th*:  entire  tro|iosphere  at  the  s(iecified  elevation  angle.  Then,  If  the  number  punchc'd  in 
Column  74  of  the  numerical  data  in|)ut  card  was  5  or  less,  the  calculated  range  or  ranges 
are  printed  for  the  Swerling  case  or  cases  specified.  Along  wi'li  each  range  figure  are  also 
given  the  calculated  signal-to-noise  ratio  in  decibels  and  the  trcjiospheric  attenuation  for 
that  range. 

If  the  “cafe”  jiaramc'ter  of  Column  74  is  6  or  7,  the  printed  output  is  modified  slightly 
to  refiect  the  fact  that  the  range  has  been  calculated  on  the  basis  of  an  assumed  signal-to- 
iioise  ratio  rather  than  for  specific  probabilities  of  detection  and  false  alarm  and  a  specific 
fluctuation  model. 

Figures  1  through  6  are  illustration.s  of  the  input  data  cards  and  resultant  printed 
output  for  three  different  'case"  options,  namely  1,  5,  and  7.  Cases  0,  2,  3,  and  4  produce 
output  results  similar  to  that  shown  for  Case  1,  and  Case  6  produces  an  output  similar  to 
that  of  Cast  7,  except  that  the  range  is  then  calculated  for  only  one  signal-to-noise  ratio. 
The  radar  parameters  of  these  sample  calculations  are  fictitious. 


EQUATIONS  AND  ALGORITHMS 


It  has  been  mentioned  that  E<).  (12)  of  NRL  Report  6930  (1)  is  the  basis  of  Program 
RGCALC.  The  equation  is 


«mas  =  129-2 


1/4 


(3) 


The  symbols  in  this  equation  have  been  previou.sly  defined  in  this  re|)orl  except  for  those 
that  follow: 

Fi,  F^  —  pattern-propagation  factors  for  the  radar  transmilter-to-target  and  target-to- 
receiver  paths,  respectively.  In  Program  RGCALC,  f’,  =  F,  =  1. 

Tg  —  the  receiving  system  noise  temperature,  kelvins 


V'„  —  the  visibility  taclor,  or  pn'detection  signal -to-noise  ratio  required  for  the 
specified  [>robability  of  detection  of  the  target  echo 


NRL  REPORT  7448 


11 


L  —  the  system  loss  factor,  equal  to  the  product  of  all  the  component  loss  factors.* 
(In  decibels,  I.  is  the  sum  of  the  <-omponent  decibel  losses.) 

Program  RGC.ALC  is  an  executive  routine  which  reads  the  data  cartls  and  then  trans 
mits  the  data  to  a  subroutine  named  RANGE.  This  subroutine  calls  other  subroutines 
which  calculate  the  system  noi^c  temperature  7,,  the  signal-to-noise  ratio  V^,  and  a  table 
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Fi^.  ]  ~  Data  cards  for  calculating  ran^c  of  a  UHF  search  radar,  SwiriinR  case  1  tarKcl  fluctuation 

(1  punch  in  Col.  74  of  second  card). 


*Lo66  factor  is  sometimes  expressed  as  a  number  le.ss  than  one,  in  which  C4Se  it  would  belong  in  the 
numerator  of  the  range  equati6n.  Mere  th*-  coireci  engineering  definition  of  loss  as  the  reciprocal  of 
gain  is  employed,  so  that  L  •  I 
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RADAR  NAWE  P)R  DESCRIPTI§N 
MCTlTIOtlS  UHE  search  RADAR 


RADAR  AND  target  PARAMETERS  (INPUTS)  •• 


•  I 


PULSE  P0EER,  KW  . . 

PULSE  length,  H1CR0SEC  . . . 

transmit  antenna  gain,  Oi 

RECEIVE  antenna  GAIN,  DB 

fREOUENCY,  mhZ  . . . 

RECEIVER  N0ISE  rACTSR  (TIQURE),  D» 

pandnidth  CRRRECTieN  pactsr,  ne 

ANTENNA  eHMtC  L6SS|  DB  7  . 

transmit  TRANSMISSIBN  line  LBSSi  D8  . 

receive  TRANSMISSIBN  LINE  LBSS#  OB  7 . 

SCANNINGvANTENNA  pattern  Less,  OB  . . 

MISCELLANEeuS  LOSS,  DB  . . . 

NUMBER  er  PULSES  INTEGRATED  7  ...  . . 

probability  eP  oetectibn  ,,.,1. i,.,... 
FALSE-ALARM  PRBBABILITY,  negative  PBrER  bp  ten 
target  CR0SS  SECTI8N,  SQUARE  MfcYERS 

target  ELEVaTIBN  angle.  DEGREES  . . 

average  SGLAR  and  QALACTI-C  NeiSE  ASSUMED 
PATTERN-PR0PAGATI0N  ACT0RS  ASSUMED  *  X. 


900,  P 
20.0000 
2?.P 
22. P 
315, P 

2,5 

0,5 

ft ,  ft 
0,5 

0,5 
1,0 
0 . 0 
45 
0,500 
B,0 
5,090P 

Kso 


/ 


CALCULATED  quantities  (BUTPUTS)  •• 

NOISE  temperatures,  degrees  KELVIN 

antenna  <TA)  . . . 7...y, 

RECEIVING  transmission  LINE  CTR)/, 

receiver  (TEJ  . . 

Tb  X  LlNE*Le5S  FACTBR  *  TEI 
SYSTEM  (TA  •  TR  *  TEI)  ,  ,  ,  yO . . 

Two-way  attenuation  through  enti»e  troposphere 


DB 


251.0 

35,4 
225,  T 

253.2 

540.2 
1,0 


SWERLING 

fluctuatibn 

case 


1 


signal* 

Ti-N0ISF 

RATIO,  08  / 

✓ 


1,00 


TRePBSPHERIC 

attenuation, 

OectoELS 


O.PT 


range, 

NAUTICAL 

MILES 


209,9 


Fig.  2  —  IVogTsm  RGCALC  output  fur  daU  card,  of  Fig.  1 
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Fig.  3  —  DaU  card*  for  calculating  range  of  a  microwave  volume^scanning  radar,  Swerling 
caies  0,  1,  2,  3.  4  fluctuation  (5  punch  in  Col.  74  of  aecond  card) 


(array)  of  tropospheric  abs  .'rption  losses  in  decibels  at  range  increments  along  the  (refracted) 
ray  path  for  the  specified  elevation  angle.  Then  a  system  loss  factor  L,  with  tropospheric 
absorption  loss  omitted,  is  calculated.  The  range  equation  arithmetic  is  then  performed, 
using  the  input  data  and  the  calculated  values  of  T,,  V„,  and  L.  Then,  in  the  table  of 
calculated  absorption  loss  values,  a  value  of  this  loss  corresponding  to  the  calculated  range 
is  found  by  interpolation.  The  range  is  correctefl  by  a  factor  corresponding  to  this  loss 
factor;  then  the  new  loss  factor  corresponding  to  this  corrected  range  is  found,  again  by 
interpolation;  this  iteration  is  repeated  until  the  last  correction  corresponds  to  less  than 
0.1  dB.  in  a  subroutine  named  ITEF  AT. 
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radar  NA>iE  eR  DESCRtHTlfiN  •« 
riCTlTieuS  MIC^SwAVE  V0UL'He*SCAN  RADAR 

RADAR  ARC  Target  parameters  (inputs) 


RULSE  POnER,  Km  . . iti.i.ofT, . «...  12CQfC 

PULSE  LEnCTh,  RiCRdSEC  V .........  .  60.0000 

transmit  antenna  Gain,  OB  « 7. . . .  34, r 

RECEIVE  antenna  GAIN.  OB  .,,,.....,.7 .  34,0 

(REOuENCYi  phT  . . I . 7..., .  24S>0,C 

receiver  n91ee  TAcreH  (eioureii  os  3.0 

sANDeIOTn  CSRRECTtSN  PaCT^R,  DB  .  1.0 

antenna  RmmIC  '.SSSi  or  7  •  s,n 

transmit  TRanSmISSICn  line  L«SS,  cb  .  o,o 

RECEIVE  TRanSHISSISN  l|nE  LSSS,  DB  1,2 

SCANNiNGtANTENNA  PATTERN  IBSS.  DB  .  3,2 

MlSCEwLANEBUS  LOSS,  DR  V  .  O.T 

NUMBER  Qf  PUlSE5  INTEGRATED . .  2 

PR0BA8ILITY  er  OETECIieN  . .  O.ROO 

fA^SE-ALARM  PReBABlLlTY.  NEGATIVE  PEJmER  BP  TEN  6,0 

:aRGeT  CRQSS  SECTI9K,  SQUARE  **ETERS  ,..., .  1,0000 

target  ElEVATIQN  angle,  degrees  0,40 

MAXIMUM  SOLAR  AND  GALACTIC  NStSE  ASSUMED 
PaTTERN«PR0PaGATIPn  fiCTBRS  ASSUMED  .  1. 


calculated  quantities  (OUTPUTS)  •• 


NOISE  TEMPERATURES,  DEGREES  KELVIN  t* 

antenna  (TA)  . . . . 7 .  211,1 

RECEIVING  transmission  LINE  (TR)  92,3 

receiver  (TE)  . ,..,....,.7 .  280,6 

Tb  X  line-loss  PACTBH  a  TEl  . .  380,9 

system  (TA  •  TR  *  TED  . .  683,9 

Tue-aAY  attenuation  TWReuGM  entire  TRRPeSPMERE,  DR  3,6 


ShERLING  signal-  TRBPeSPMERlC  range, 

PLUCTuATIBN  T8-N0tSE  ATTENUATltN,  NAUTICAL 

CASE  RATIO,  OB  decibels  PILES 


0  10, 62  2,88  126.9 

1  IR,6V  2, IT  83,2 

2  14,83  2,22  101,9 

3  14,83  2.92  101,9 

4  12,66  2.69  113,0 


KiK,  I  -  Program  KUCALC  ouipui  lor  rt*ia  ca^os  oi  Kig.  .» 
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RaOAH  name  or  OESCRIHTieN  •<! 
fictitious  mJlLIRETER.WAVE  tracking  radar 

RADAR  AND  TARGET  PARAMETERS  (INPUTS)  -• 

Pulse  pomer,  kw  . . so.o 

PULSE  LENGTW,  MlCRySEP  7  ......... .  o,29on 

TRANSMIT  ANTENNA  GAIN,  0«  ,.,,,,..,.7..., .  54. fl 

RECEIVE  ANTENNA  GA|N,  06  V  ...  . .  54.0 

FRE3W6NCY,  mm2  ...,.,....7..., .  SQOOO.'^ 

RECEIVER  NeiSE  PACTOH  (PIGUREJ,  D8  . .  15,0 

BANDmIDTM  C0RPECTI8N  PACTflH,  DB  0.9 

ANTENNA  OWMIC  |,6SS,  08  7  ...  . .  0.1 

transmit  transmission  line  Less,  DB  0,y 

RECEIVE  transmission  LINE  L0SS,  Ob  V . . 

SCANNING«ANTENNA  pattern  L6SS,  OB  0.9 

MISCELLANEOUS  L^SS,  DR  7 .  0.3 

SIGNAL-TO^NPISE  RATje,  DB  V  ...  . .  3,0 

SIGNAL-TB-NQISE  RATI9,  DB  ,..,,.,.,.7...,. .  0,0 

target  cross  SECTION.  SQUARE  METERS  . .  0.2000 

target  elevation  angle,  degrees  . .  0,00 

MINIMUM  solar  and  GALACTIC  NOISE  ASSUMED 
PATTERN-PROPAGATIRN  FACTORS  ASSUMED  ■  I, 


calculated  quantities  (OUTPUTS)  •• 

NOISE  temperatures,  DEGREES  KELVIN  r, 

antenna  (TA)  . ,..,.,..,.7,..,..,...  ?S2,4 

RECEIVING  transmission  LINE  (TR)  . .  191,3 

receiver  (TE)  . .  8080,6 

TE  X  LINE-L8SS  FACTOR  ■  TE  I  . . .  14738,1 

SYSTEM  (TA  *  TR  *  TED  . .  15212,0 

TW0-hAY  ATTENUATION  THROUGH  ENTIRE  TROPOSPHERE,  D8  34,4 


RANGE  *  22,1  N,  MI,,  TRSPeSPHERJC  ATTENUATIBN  •  7,27  DB 

FOR  SPECIFIED  SIQNAL’Te-NOISE  RATie  f  3.00  DB 

range  *  25,0  N,  Ml,,  TRSPeSPHERIC  ATTENJATieN  •  8,15  DB 

FOR  SPECIFIED  S I  ON AL -TO-NO I SE  RATI9  f  O.OO  DB 

Fig.  6  —  Piogram  RGCALC  output  for  data  card,  of  Fig.  5 
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Figure  7  is  a  flow  chart  showing  the  sequence  of  events  in  the  computation.  As 
indicated,  checks  are  made  at  several  points  whicli  permit  certain  portions  of  the  calculation 
to  be  omitted  if  more  than  one  range  calculation  is  being  made  and  if  results  of  some  steps 
in  the  preceding  calculation  can  be  used. 

In  the  fol’owir.o  sections,  the  algorithms  of  the  various  subroutines  of  the  program 
will  be  described.  The  names  in  parentheses  following  the  subroutine  name  are  its  parameters 
in  the  culling  sequence. 


Subroutine  ALPHA  (FMHZ) 

This  subroutine  performs  the  initial  step  in  the  computation  of  tropospheric  absorption 
loss  and  noise  temperature;  it  computes  a  set  of  absorption  coefficients  in  decibels  per 
nautical  mile  for  a  set  of  attitudes  above  sea  level  from  0  to  100,000  ft.  The  ‘‘irst  21  of 
these  altitudes  (from  0  to  2000  ft)  are  at  intervals  of  100  ft;  the  next  28,  to  30,000  ft, 
are  at  intervals  of  1000  ft;  the  next  20,  to  70,000  ft,  are  at  intervals  of  2000  ft;  and  the 
last  6,  to  100,000  ft,  are  at  intervals  of  5000  ft.  This  graduation  of  height  intervals 
reflects  the  fact  that  the  absorption  coefficient  changes  more  rapidly  in  the  lower  atmosphere 
than  it  does  at  higher  altitudes. 

The  absorption  coefficients  are  calculated  at  each  of  these  altitudes  fo*'  both  oxygen 
and  water  vapor,  and  the  two  coefficients  are  added  to  obtain  a  total  absorption  coefficient. 
The  resulting  array  of  75  coefficients,  for  the  frequency  FMHZ  (first  parameter  of  the  call¬ 
ing  sequence)  is  named  ALPH  tJ,  75),  J  =  3.  It  is  transmitted  as  output  via  a  COMMON 
block  named  PTR. 

The  computation  is  done  using  the  theory  of  Van  Vleck  as  described  in  NRL  Report 
7010  (1),  except  that  in  the  previous  version  of  the  subroutine  described  there.  Van  Vleck ’s 
centroid  approximation  was  used,  and  it  is  not  valid  in  the  region  near  the  oxygen  resonances 
from  about  50  to  70  GHz.  The  new  version  of  the  subroutine,  as  now  used  in  Program 
RGCALC,  performs  a  more  exact  calculation  by  summing  the  separate  absorption  contribu¬ 
tions  of  each  of  46  individual  oxygen  resonance  frequencies.  Consequently,  range  calcula¬ 
tions  can  now  be  made  within  as  well  as  outside  the  frequency  region  50  to  70  GHz. 

The  calculations  are  made  for  the  standard  dry  atmosphere  known  as  the  U.S.  Exten¬ 
sion  to  the  Internationcl  Civil  Aviation  Organization  (ICAO)  Standard  Atmosphere  (10). 

The  model  tor  water-vapor  content  of  the  atmosphere  is  based  on  a  humidity  profile 
given  by  Sissenwine  and  others  (11)  as  representative  of  the  midlatitude  mean  humidity. 

This  midlatitude  mean,  however,  has  a  surface  water-vapor  content  of  5.947  gm/m^,  and 
it  was  desired  (in  accordance  with  prevailing  practice)  to  adopt  the  surface  water-vapor 
density  value  of  7.5  gm/m^  for  the  absorption  computations.  Therefore,  the  values  of 
the  Sissenwine  model  were  all  multiplied  by  the  factor  7.5/5.947  =  1.261.  The  tabulation 
given  in  Sissenwine’s  report  is  for  altitudes  at  intervals  of  2  krn  in  the  region  of  interest 
here.  The  values  corresponding  to  the  altitudes  specified  in  Subroutine  ALPHA  were 
obtained  by  means  of  en  interpolation  technique  developed  for  digital-computer  plotting 
of  a  smooth  curve  through  a  set  of  data  points  (12). 

The  sets  of  75  values  of  pressure,  temperature,  and  water-vapor  density  values  defining 
this  model  atmosphere  me  eritered  into  Subroutine  .''.LPH.A  in  the  form  of  Fortran  DATA 
statement*:  (arrays  PP,  TT,  and  RR),  thus  obviating  any  necessity  of  reading  them  in  from 


Flow  chart  for  Program  RGCALC 
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a  separate  deck  of  data  cards  each  time  the  program  is  run.  The  set  of  46  oxygen  reso¬ 
nance  frequencies  is  similarly  entered  via  DATA  statements;  these  frequencies  are  separated 
into  two  arrays,  FTRP  and  FTRM,  cuiTesponding  to  two  classes  of  quantum-mechanical 
state  transitions  of  the  oxygen  molecule.  The  details  of  the  calculations  are  described  by 
Meeks  ano  Lilley  (13);  their  formulations  of  the  oxygen  absorption  equations  were  employed, 
(The  frequencies  FTRP  correspond  to  their  symbol  e,v*>  FTRM  corresponds  to 

The  only  deviation  from  the  Meeks  and  Lilley  calculations  was  the  use  of  a  slightly 
different  dependence  of  line  width  on  altitude.  The  model  of  Reber,  Mitchell,  and 
Carter  (14)  was  used  for  this  pail  of  the  calculation. 


Subroutine  ATLOSS  (FMHZ,  ELEV,  ATMP) 

The  input  parameters  are  FMHZ,  frequency  in  megahertz,  and  ELEV,  elevation  angle 
in  degrees.  The  output  parameter  ATMP  is  the  tropospheric  noise  temperature  computed 
for  the  specified  frequency  and  elevation  angle.  Another  output,  transmitted  via  COMMON 
block  RGA,  is  an  array  of  absorption  values  named  ATTN  (decibels)  corresponding  to  a 
set  of  range  values  along  the  ray  path  at  angle  ELEV,  corresponding  to  the  altitude  values 
of  Subroutine  ALPHA.  The  corresponding  array  of  range  values  RG  is  similarly  transmitted 
as  output. 

Subroutine  ALPHA  is  called  by  Subroutine  ATLOSS,  and  the  resulting  airay  of  75 
values  of  absorption  coefficients  ALPH(J,75),  J  =  3,  is  used  to  calculate  cumulative 
absorption  along  the  ray  path  at  angle  ELEV.  The  ray  path  in  the  refracting  atmosphere 
is  computed  by  numerical-integration  ray  tracing,  from  the  formula  (15) 


(4) 

in  which 

dR  _  nf-h) 

,/~  [noCosOola 

[n(l.h/r„)J 

(5) 

Here  R  is  the  radar  range  corresponding  to  height  /ij  as  measured  along  the  ray  path  of 
the  initial  elevation  angle  ,  n(.'i)  is  the  refractive  index  height  profile,  is  the  value  of 
n  at  /]  =  0,  and  is  the  radius  of  the  earth  (more  specifically,  it  is  the  distance  from  the 
earth’s  center  to  the  initial  point  of  the  ray). 


The  attenuation  is  then  computed  along  this  ray  path  by  numerically  integrating 
the  equation 
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in  which  =  R{h^),  and  -){h)  is  the  absorption  coefficient  at  height  h,  as  given  by  ihc 
array  ALPH  (J,  76),  J  =  3,  The  derivative  ds/dh  is  equal  to  [l/n(/t)l  •  dRjdh',  that  is, 

R  is  the  radar  range  measured  along  the  ray  path,  and  s  is  the  geometric  distance  along 
the  same  path.  The  derivatives  dsjdh  and  dRidh  are  computed  in  a  short  subroutine 
named  DDH{H),  in  which  the  single  parameter  is  the  height  H.  The  derivatives  are  named 
DSDH  and  DRDH  and  are  transmittf^  to  Subroutine  ATLOSS  via  a  COMMON  block 
named  DRS. 


The  refractive  index  model  used  is  given  by  (2) 

n(/j)  -  1  >  0.000313  e-''''  (7) 

where  k  =  4.3848  X  10~®  if  h  is  in  feet.  The  earth’s  radius  is  assumed  to  be  6370  km,  or 
2.0899  X  10'^  It. 

A  special  technique  is  used  to  perform  the  integrations  in  the  vicinity  of  h  =  0  for 
the  special  case  0^  =  0,  because  in  that  case  dR/dh  and  dsjdh  both  become  infinite.  This 
technique  was  described  by  the  author  in  a  paper  published  in  1968  (15). 


The  tropospheric  noise  temperature  (Fortran  name  ATMP)  is  computed  by 
numerical  integration  of  the  equation 


‘  atm 


0.2303  i(R)  Tf(R)  exp 


-0.2303  7(r)c/r 


dR. 


(8) 


where  dR  is  taken  alor®  the  ray  path.  T,  is  the  thermal  temperature  of  the  troposphere; 
its  values  are  transmitted  to  Subroutine  ATLOSt  from  ALPHA  via  COMM(JN  block  TMP. 

The  previously  described  modification  of  Simpson’s  rule  cannot  be  used  to  perform 
this  integration  because  it  is  an  integration  with  respect  to  R  (range)  rather  than  h  (height). 
The  h  intervals,  as  described  in  the  section  on  Subroutine  ALPHA,  are  uniform  (over  each 
of  the  four  height  regions);  however,  the  corresponding  R  intervals  are  not  uniform. 
Another  special  modification  of  Simpson’s  rule  was  devised  to  handle  this  problem;  it  is 
embodied  in  Subroutine  INTGRT,  which  is  called  by  ATL(5SS  to  perform  the  numerical 
integration  of  Eq.  (8).  Further  details  of  the  absorption  and  noise  temperature  calcula¬ 
tions  will  be  given  in  a  report  to  be  written  in  the  near  future,  in  which  curves  for  absorp¬ 
tion  and  noise  temperature  as  functions  of  frequency  and  elevation  angle  will  be  presented. 

Subroutine  TEMP  (FMHZ,  ELEV,  K,  ANF,  ALA,  ALR,  ATMP, 

TA,  TR,TE,TEI,TSYS) 

The  input  parameters  FMHZ,  ELEV,  and  ATMP  are  the  same  as  those  of  Subroutine 
ATLOSS.  The  other  input  parameters  are 

K  —  galactic  noise  code  (—1,  C,  -i-l) 

ANF  -  receiver  noise  factor  NF,  dB 
ALA  "  antenna  loss  factor  L^,  dB 
ALR  —  receiving  line  loss  factor  L^,  dB. 
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The  output  pnrameters  are 

TA  —  antenna  noise  temperature  ,  kelvins 

TR  —  receiving-transmission-line  noise  temperature  7^ 

TF  —  receiver  noise  temperature  Tf. 

TEI  —  product  of  Tp  and  L, 

TSYS  —  system  noise  temperature  Tj.  , 

The  antenna  noise  temperature  is  computed  by  use  of  equations  given  in  NRL  Report 
6930,  p.  49,  and  Report  7010,  pp.  40  through  44  (1).  The  sky  temperature,  named  TAl, 
is  first  computed  from  the  equation 

Tsky  =  (Tgsl  *  7'66)//>atm  +  ^a(sun)  (9) 

in  which  Tggi  is  the  galactic  noise  temperature,  T{,{,  is  the  cosmic  blackbody  temperature 
(2.7  K),  is  tl  -■  atmospheric  loss  factor  (expressed  as  a  power  ratio  >  1),  Tg(sun) 

the  solar  contribution  to  the  antenna  temperature  (assuming  that  the  sun  is  in  an  average- 
level  sidelobe  of  the  antenna  pattern),  and  Tatm  is  the  atmospheric  noise  temperature 
(ATMP,  obtained  from  Subroutine  ATLOSS).  The  galactic  temperature  is  given  by 

7’gal  'T'lOO  •(100//^MHz)^  ‘^-  (10) 

The  quantity  Tjog  is  the  galactic  temperature  at  the  frequency  /'mh/,  “  100  MHz.  Its 
numerical  value  depends  on  the  galactic  noise  code  K  according  to  the  following  prescription. 


K 

Tjoo  (kelvins) 

-1 

500 

0 

3050 

+l 

18,650 

The  solar  contribution  to  antenna  temperature  is  obtained  from  a  table  of 

values  of  the  solar  noise  temperature  entered  via  a  DAT.A  statement;  the  table  corresponds 
to  frequencies  in  the  range  100  MHz  to  10  GHz,  and  the  values  are  taken  from  Fig.  6  of 
Report  7010,  p.  43  (1).  At  frequencies  between  the  tabulated  values,  Tgu,,  is  found  by 
linear  interpolation.  Above  10  GHz,  Tjun  'S  assumed  to  have  the  constant  value  10,000 
kelvins.  The  solar  temperature  Tgun 's  related  to  Ta(sun)  hy  the  equation 

7'a(sun)=7’gun><  4.7b  X  (IT) 

The  numerical  factor  takes  into  account  the  assumed  unity-gain  average  sidelobe  level,  the 
ratio  of  the  sun’s  noise  diameter  to  the  total  solid  angle  {47i  steradians)  viewed  by  the  antenna 
including  its  side  and  back  lobes,  and  the  assumption  that  the  sun  is  on  the  average  ten  times 
noisier  than  indicated  by  the  referenced  curve,  which  portrays  the  “quiet  sun.”  Then, 
is  decreased  by  a  factor  of  10  if  K  =  —1  and  increased  by  10  if  K  =  +1 ,  where  K  is  the  galactic 
(and  solar)  noise  code. 
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The  resulting  value  of  sky  temperature  TAl  is  multiplied  by  0.876  to  take  into  account 
the  fraction  of  the  total  antenna  pattern  subtended  by  the  sky,  and  to  this  is  added  the 
contribution  due  to  antenna-loss  noise,  in  accordance  with  Eq.  (37)  of  NRL  Report  6930. 

The  transmission  line  and  receiver  noise  temperatures  Tr  and  Te  are  computed  in 
accordance  with  Eqs.  (40)  and  (41)  of  NRL  Report  6930,  n.  50,  and  combined  to  give 
the  system  noise  temperature  by  the  equation 

T,  =  7,  (12) 

where  is  the  receiving  line  loss  factor.  The  product  is  also  reported  to  Subroutine 
RANGE  as  the  parameter  TEI,  and  is  printed  as  an  intermediate  output  of  the  program 
along  with  T^,  T^,  T^,  and  T,. 

Subroutine  PDSN  (PDT,  PEA.  NPULS,  KASE,  SDB) 

The  purpose  of  this  subroutine  is  to  find  the  signal-to-noise  ratio  required  for  detection 
SDB,  for  the  specified  probability  of  detection  PDT,  false-alarm  probability  PFA  (expressed 
as  a  positive  number  representing  the  negative  power  of  ten),  number  of  pulses  integrated 
NPULS,  and  Swerling  fluctuation  case,  KASE.  Subroutine  PDSN  does  not  perform  the 
calculation;  it  merely  manages  it  by  calling  other  subroutines.  The  actual  calculation 
requires  an  iteiation,  v, liicii  it  performed  by  Subroutine  INVERS,  called  by  PDSN.  Before 
calling  INVERS,  PDSN  estimates  a  range  of  decibel  values  (lower  value  DBl,  uppei  value 
DB2)  likely  to  contain  the  true  value  SDB.  An  empirical  formula  is  used  for  this  piupose. 
This  procedure  minimizes  the  number  of  iterations  required.  Subroutine  PDSN  is  called 
from  Subroutine  RANGE,  and  the  resultant  value  of  SDB  is  used  as  a  factor  in  the  range 
calculation.  When  the  “case”  parameter  of  Subroutine  RANGE  is  6  or  7  (Col.  74  of 
the  data  card),  PDSN  is  not  called,  since  the  signal-to-noise  ratio  is  then  a  direct  input 
and  need  not  be  calculated. 


Subroutine  INVERS  (XMIN,  XMAX,  XL(9,  XHI,  NSIG,  LIM, 

NCd,  X,  FI,  FT,  F) 

This  subroutine  performs  an  iteration  to  determine  the  value  of  the  argument  X  of 
a  function  F(X)  which  will,  within  a  specified  accuracy,  cause  F(X)  to  equal  FT,  a  specified 
value  of  the  function.  The  accuracy  parameter  is  NSIG  —  the  number  of  significant  figures 
to  which  agreement  is  desired  between  F(X)  and  FT.  LIM  specifies  a  limit  on  the  number 
of  iterations  fjermitted,  and  NOI  (output  parameter)  reports  the  number  of  iterations 
actually  performed.  FI  is  the  actual  final  value  of  F.  It  is  required  that  F  be  a  monotonic 
function  of  X  within  the  permissible  range  of  variation  of  X,  which  is  from  XMIN  to 
XMAX,  and  that  tne  value  FT  exists  within  this  range.  The  input  parameters  XL(^  and 
XHI  define  a  ’•egion  in  which  it  is  guessed  that  the  desired  value  of  X  will  be  found.  If 
no  knowledge  exists  by  which  to  e.stimate  this  region,  XL©  and  XHI  can  be  set  equal  to 
XMIN  and  XMAX;  however,  the  mere  narrowly  the  region  is  defined,  the  fewer  will  be 
the  iterations  required. 

If  the  slope  dF/'dX  becomes  nearly  zero  in  some  part  of  the  range  from  XMIN  to 
XM.\X,  and  if  this  region  is  contained  between  XL©  and  XHI,  the  iteration  will  converge 
very  slowly,  or  conceivably  not  at  all.  It  is  for  this  reason,  as  mentioned  earlier,  that 
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values  of  probability  of  detection  greater  than  0.99  should  not  be  specified.  Evidently 
Subroutine  MARSWR,  which  is  involved  in  the  iteration,  does  not  define  well  the  slope 
of  the  “function”  in  this  region. 


Function  PD  (SNDB) 

In  order  to  define  a  function  on  which  INV'ERS  can  operate,  the  Fortran  FUNCTION 
PD  is  used,  with  signal-to-noise  ratio  (dB)  as  the  argument  SNDB.  This  function  merely 
calls  Subroutine  M.ARSWR,  which  calculates  the  probability  of  detection. 


Subroutine  MARSWR  (SNDB,  N,  FA,  KASE,  PN) 

As  has  been  discussed,  this  subroutine  is  basically  the  subroutine  of  Fehlner  (6)  which 
he  mmod  MARCUM.  It  was  renamed  MARSWR  after  a  few  changes  in  it  were  made  to 
adapt  it  to  the  requirements  of  Program  RGCALC.  The  input  parameters  are  SNDB, 
the  signal-to-noise  ratio  in  decibels;  N,  the  number  of  pulses  integrated;  FA,  the  false-alarm 
exponent;  and  KASE,  the  Swerling  fluctuation  case. 

As  has  been  mentioned,  Fehlner’s  original  subroutine  calculates  the  probability  of 
detection  PN  for  a  specified  value  of  Marcum’s  false-alarm  number  rather  than  on  the 
basis  of  false-alarm  probability.  The  power-of-ten  exponent  of  the  false-alarm  number  is 
named  FAN  in  the  subroutine.  The  relation  between  FAN  and  the  false-alarm  probability 
Pfa  is 


FAN  =  logio 


logf  0.5 
>ogc(l  -Pfa) 


(13) 


This  relationship  is  used  in  Subroutine  MARSWR  to  convert  the  input  parameter 
FA  to  the  internal  parameter  FAN. 

In  using  Subroutine  MARCUM,  it  was  found  that  an  appreciable  portion  of  the 
computing  time  is  spent  in  computing  the  bias  level  YB.  If  succesive  calls  to  MARCUM 
are  made  with  the  same  values  of  N  and  FA  (but  with  different  values  of  SNDB  and  K.ASE), 
it  is  not  necessary  to  repeat  the  bias-level  calculation.  Therefore  a  provision  for  omitting 
that  part  of  the  calculation,  when  succcs.sive  calls  to  MARSWR  are  made  with  the  same 
values  of  N  and  FA,  has  been  added  to  the  subroutine. 

Functions  named  DGAM,  DEVAL,  GAM,  and  SUMLOG  are  part  of  the  M.\RCUM 
subroutine,  which  was  originally  written  in  Fortran  II  for  use  with  an  IBM  computer. 

These  functions  are  also  incorporated  into  MARSWR.  The  only  changes  made  in  them 
were  those  necessary  to  adapt  them  for  use  on  the  NRL  CDC-3800  computer.  (Some  of 
these  adapting  changes  were  made  by  Stanley  Gontarek,  of  the  Naval  Air  Systems 
Command.)  A  further  slight  change  was  made  in  Function  SUMLOG;  the  array  named 
A  therein  was  given  a  dimension  1000  rather  than  200  as  in  the  original  MARCUM  sub¬ 
routine.  This  increased  dimensioning  saves  computing  time  if  successive  calls  to  SUMLOG 
are  made  with  values  of  N  greater  than  200. 


24 


L  V  BLAKE 


FORTRAN  PROGRAM  LISTING 

The  Fortran  program,  subroutines,  and  functions  are  listed  on  the  following  pages. 

The  names  of  the  listed  routines  and  their  computer  lengths  (numl)er  of  locations  required) 
are  as  follows: 


Name 


Program  RGCALC 
Subroutine  RANGE 
Subroutine  TEMP 
Subroutine  ATLiZiSS 
Subroutine  DDH 
Subroutine  ALPHA 
Subroutine  INTGRT 
Subroutine  ITERAT 
Subroutine  PDSN 
Function  PD 
Subroutine  INVERS 
Subroutine  MARSWR 
Function  DGAM 
Function  DZVAL 
Function  GAM 
Function  EVAL 
Function  SUMLOG 


Locations  Required 

Octal  Decimal 

232  154 

2030  1048 

362  242 

756  494 

112  74 

520  336 

211  137 

220  144 

230  152 

63  51 

714  460 

3123  1619 

216  142 

lOG  70 

206  134 

111  73 

4203  2179 

16525  7509 


Totals 
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SuPReuT  INE  R*NGb(PT«<h,TAUNS|OTOB,CROe,S10S'*.»^MW2i*L*PBi  ALWf’Bi 
1  »lTI;0;  *uPDR,  *t*CR,Cdn8.*Nf  08,N*ULS,Pn,Pr*,'‘*SE,El.EV»*Jef  $«  ) 

CeMMQN/RGA/  RG(T5).  *TTN<3,79) 

n  (kase  ,eo.  5)  in,ii 

10  IkaSE  ■  0 

ce  TR  12 

11  IkaSE  ■  KASE 

12  call  lEMPCf’MM2»ELEV,NttlSEiANr0R,AL*06,4LRDP.*TMP,Ti,TR,TE,T6I, 
1  TSYS) 

n  (KASE  ,GE.  6)  60,61 

60  SNDBtPO 
GS  TR  62 

61  call  PnSNCPD.PfA.NPULSi IKASEjbSOal 

67  r ACDa»GTOa*GRUR»Cbr>a« ALTnu^ALPObaALXDBuSNDR 

PACilO.*»(P*CDP»,l) 

SNr,t]«i29 .2*  <PTKh«TA|JRb«S16SR*E  IC/  (rMh2*rMM2*TStS)  ,25 
RNGAtRNGO 

C*I.L  |TEfiAT(BNCA,ATT) 


PRINT 

1 

print 

lOOi 

PTKW 

print 

lOli 

Tau-S 

print 

102i 

fiTDP 

print 

103, 

CRoe 

print 

10«i 

P-hZ 

PRINT 

10^1 

ANfOB 

print 

106, 

CROP 

print 

10»i 

AlaOB 

PRINT 

106| 

AlTOB 

print 

110, 

AlROB 

print 

nil 

ALPDB 

print 

112| 

AlXDB 

(KASE  , 

CE.  6| 

63,6 

63 

PRINT 

ll7, 

Snob 

IF  (KASE  , 

EO.  71 

PRIN 

ce  TO 

65 

64 

print 

107, 

NPLILS 

print 

Il3, 

PO 

print 

116, 

pr  A 

65 

print 

115, 

SIGSN 

print 

116, 

ElEV 

IF  (NelSEI 

50,51, 

52 

50 

print 

55 

Ge  re 

53 

51 

print 

56 

ce  re 

53 

52 

PRINT 

57 

53 

print 

56 

print 

60 

print 

2 

PRINT 

120 

print 

121, 

TA 

print 

122, 

TR 

print 

123, 

TE 

PRINT 

1123 

.  TEI 

print 

126, 

TSYS 

print 

125, 

ATTNI3 

,75) 

26 


L.  V.  BLAKE 


05/12/72 


PROBNAM  R6C*LC 

OlMCf^SION  NAMC(IO) 

2  READ  3*  name 

ir  (EOF«60)  lOtll 
10  STOP 
n  RSINT  22 

22  P0RmaT<15«»»?'*0AR  NAME  OR  OCSCRlRTlOw  -•  •/» 

PRINT  30*  NAME 

READ  A*  PT*TAU*0T*6R*S10«RN*AtA*ALR»ALT*ALRtALX*ce«ANF,NP,P0* 
1  fa»ka*el*ns 

CALL  Range (PTiTAU*OTtGR«SlO*rM,ALA*ALR» alt. ALR»AL*»CB*ANF,NPtPO* 
1  P*tHA«EL*NS) 

PRINT  5 

no  TC  2 

1  format («X«I4) 

3  format (lOAB) 

30  FORMAT ( lOXtIOAG//) 

4  format (2F«.0«2F4.0«2F6.0t7FA*0»lS«2F4,0» 11 fFA.O, 12) 

5  format <IM1) 

FNO 
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\f  {KASE  ,0E,  6)  *A,«/ 

•«  PRINT  1««|RN0A,ATT 
PRINT  IAT,  snob 
IF  (KASE  ,E0.  TJ  JTO.  m 
1^0  SNOtl  •  PFA 

DIFF  f  SNPevSNOHl 

PNOI  •  PNSO«10.P«(D|Ff *,029) 

caul  |TEPAT<RnG|.ATT) 

PRIN^  iAA|  RNai,  ATT 
PRINT  tA7|  SN0B| 

IM  RCTyRh 
AT  PRINT  150 
print  151 
print  152 
PRINT  153 

PRINT  14,]KASE.iN09«AtT.RN0A 
IF  (K4SE  ,EQ,  5)  20,21 
2C  oe  30  I«1|A 

CALL  P0SN<PD,PFA,NPULS# 1 f SN091 1 
OlFPtSNDBaSNOBl 
FAC»10.«»(OIFI •,02B) 

RN0|«RNG0«FAC 
CALL  ITeRAT(9NO| , ATT) 
print  14.  I,SnDBIiATT.PnGI 
JO  CeNTiNUE 
21  RtTyRN 

1  reRHAT  (ISKaPADaR  and  target  parameters  IINPuTS)  •••/) 


100  FeRMAT<l5)(»PULS|  PBwER#  RW  I  . . .  Fn  ,  l ) 

Ifal  F6RMAT(15X*PULSE  L^NCTN.  HIQRBSEC  >•••>•<>**^^*>^1 

102  fbrmat<i5a»tran5mit  antenna  gain,  ne  .••fh.d 

103  FeRMAT(l9x*RECE|Vk  ANIENNA  GAIN,  OB  . Fn , 1 ) 

lOA  reR*lAT<l5x«FR60yENCT.  mmT  •»Fil,l) 

105  F8RNAT<l5x*Rec61VER  NelSE  FACTBR  (FIGgPE),  C3  . ‘.Fil,!) 

lOA  FBRMATI i5*«BanDnIOTr  CPHRbCTIBN  FACTSR,  DB  . . Fl 1 , 1 ) 

ICT  F6R«AT(l5x*NOMflfeR  RF  PULStS  INTEGRATED  . AjAx.IS) 

108  FeRMATIl5X»ANTENNA  RHMlC  LOSS.  08  . . Fn  ,  i ) 

109  feRMAT(l5x*TRANSMlT  THANSmISSION  LINE  LBSS.  CB  Fn  ,  1 ) 

110  rePMAT<l9x»RECEIVt  TRANSMISSION  LINE  L8SS,  CB  . •,Fii,i) 

111  FeRMAT(l5X«SCANN!NG*ANTENNA  PATTERN  L0SS.  D8  . F; 1  , 1 ) 

112  FePMAT<i5X»MlSCELL*NEWUS  LOSS.  D§  . Fll . 1 > 

113  F6RmAT(i5x»PR08a8ILITY  BF  OETECTIBN  .  . . . . .•.Fll, 3) 

114  FeRMAT<l1X»FALSE»AL4PM  PRBOAfllLlTY,  NEGATIVE  POwER  BF  TEN  •.Fn.l) 

115  reRM4T(i5X»TARGfcT  CRBSS  SECTI0N,  SOyARF  METERS  . . ‘.Fn.A) 


11?  FeRMAT(i9x,«SIGNAL-Tf-NeiSE  RATIB,  D8  . . . . ‘Fll,!) 


55  FeRMAl(l5X»MlNlMUM  SOLAR  AND  GALACTIC  NBISE  ASSUMED*) 

56  FSRMAT(15X*Av£RAGE  SBLAR  AND  GALACTjC  NBISE  ASSUMED*) 

5?  FeRMAT(i5X*MAXlHUM  SOLAR  AND  GALACTfC  NBISE  ASSUMED*) 


58  FBRmaT < l5X*Pi I TERN-PPBPAGATICN  FICTBRS  ASSUMED  *  1,*) 

An  . . .  /) 

2  F6RMAT  (  i5x*CAlCULATED  OL'AnTITIES  (BgTPuTS)  -.*/) 

J,i0  FePMAT(l5X*N0lSE  TEMPERATURfcS,  DEGREES  KELVIn  •-  *) 

121  rSPMATdSX*  XNTbNNA  <TA)  . .  F  1 1 , 1  ) 

122  reRMATdSX*  BECblViNC  transmission  line  (tR)  . •»Fll.l) 

123  FSRMATdSx*  RECEIVER  (Tb)  . . . •.Fll.l) 

11<3  FBRHATdSX*  TE  X  LlNE*LeSS  FJCTBP  *  TEI  r  F}  1 , 1 ) 
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*  ’'*  •  . '.rn.D 

125  r|BMiT(l5>»TW8-wAY  attenuation  TWReUCH  JNTIhE  THOboSBMERP ,  D9  *, 

1  F  7  1 1// ) 

190  rSRNAT(l5K»  SwERcl*<6  SlUNA^.  T*|RSS»>hER I C  RANGE,*) 

1»1  r0RHAT(i5X*  ri.UCTU*T!«N  T§.N6|SE  ATTENUATIBN,  nautical*) 

152  ri(*HAT(i5**  CASE  RAT18,  CS  DECIBELS  niiES*) 

193  riBNAT<l5K*  . . .  . . I. . ,,, 

14  Fe»NATI20X, ll,9X,E6,2.  9*.F6 ,  j 

16«  reBMAT(/iSX*B*NGE  •  •.fe.l.*  N,  *<1.,  TRBP9SPHER1C  aTTEN'JATIBn  .  • 

1  |F0.2<*  ue*/) 

Id?  reRMATli5X,*reR  SPECtMfcO  Sl(iNAL*T0.N*tSE  RATIO  ■  •,  F4,?,*  D8*//) 
END 


o  o  o  o 
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<f*’MtifcEV|R|ANf  tALA|ALR>ATMP,rA»Tt«Ti,TI]«TSyS) 

ttALACTiC  PttOf  *<*0  IS  fill  AVEtAOE  CiLACriC  NAtSA 
“laHEST  tEMRERATURES)!  AND  Sun  nsise 
Jr*  OUfE^  tfcVEU.  N«l  !S  r§R  MAXiNLiM  GALACTIC  NltSE 

DiifilsieN  tJJJoJ 

CiNH«N/R0A/*«g(75),  ATTN(S,79) 

01NENSI8N  rR(|),TS(e) 

OATMfRtioOl .200. #300.. 400. *900,, lOOO;, 3000. #10000, )» (TSal.ieA, 

data  JTioo  0  500.i3o9o.,l«6»0,> 

data  (TfiLKeVaJiT) 

data  (PUASTfO. ) • (ELASTaiOO. ) 

'[aJtI?5m2°'^'‘**^  eWCV.Ea.ELAST)  69  Tfi  90 

ECASTiElEv 

CALL  ATLeSS(rMM2,ELev.ATMP) 

D§  10  ta2,8 

If  (rNH2.rR(I))20,30.10 
20  JaUl 

(FR(n.rR(j))  a  TS(J> 

30  TSUN  A  TS(I) 

08  T8  40 
l«  CONTINUE 
TSUNal.OEA 

40  ATTalO,««|.ATTN(3|T9)«,09) 

>0  TASgN|4,7>E*5«TSUN*(lu,*aK)  a  ATT 

TAlaCTlOO{K*2»4(iOO,/fPHJ)aa2,9*TBLK0Y)aATT*TASUN*ATPP 

4LAAal0,aa(,ALAa,l) 

TAa  (,874aTAl  •  294,)«ALAA  a  28o. 

ALRRt;o,«a<ALRa,l) 

TRa{ALRR.i,),298, 

ANFFai0,»a{ANF«,l) 

TE»(ANFr-l, )a29o, 

TfclaTfcaALRR 

TSrs  *  TA  a  TR  a  TEI 

end 
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SUBROUT  t  NE  ATkiSS (rHH2 iElEVi A  ) 

C8HMBN  /PTB/  PP(7B),  T7(7»),  RR (75 ) , A^PH < 3 . 75 > 

CBMMON/RGA/  R6(75).ATT»4(S»75t 
DIMENSIBN  I66(4)i  DELH(4) 

CBHMON  /RRQ/RErgiRAD,QRAOiU 

CBHhBN/ORS/  0SDM3|DRDhSiAN 

DATA(RO(l)aO. ) • (ATTN(lil)aOi )« (A7TN(2i1}«0i >• ( ATTN(3i1)«0 • ) 
DATA  (REP 0S|0003i3)| (RA0a20BR5R50 tl3> I (ORAD*. 00004365) 
DATACrtASTaQ. )• tEkA!iT«iO0i)i(CBNST«,2302Sa5) 

DATA  ( t66610,x4,i0t3). (OEUHfioO, tlCOO;.2bOO. lOeOOi ) 

AT Ti(YY)aEAC2«(i. 25*71*2. •Y2*,254tV) 
ATT2(YY)*rAC2*(*.25*Vl*2.*V2*l,25*VT) 
RGi(DR)«PaCi*(1,2S*UR0Mi«2.*DRDH2>.25*0R) 
RS2(0R)*PACl*(<,25*0RUNi*2,*DRDH2*l,25*DR) 
(rHH2.E0,r|.AST,AND.PLEV,fca,EkAST)  RETURN 
EkAST  •  EkEV 
ISlOt; 

ThETA|ELEv/57, 2957795 
SN*SIN( THETA) 

CSaCeS(THETA) 

SS*SN*SN 

flPl*l,*REfO 

U»<RPJ*SN)«*2  •  2,*RE)0  .  REPOtREFO 
P  (FHMZ.EO,  FLAST)  60  TO  55 
CAkk  AkPHA(FHH2) 

FkASTlFMMZ 
55  H«0, 

RNGaO, 

ATT£Nl*ATTeN2*0, 

K*«l 

HHIN  »0, 

IF  (EkEV.EO.O.  )  MMlN*l,E-9 
call  UDH<HMiN) 

0RDHlfDRDN3 
DSt)Hl|DSDH3 
ANif AN 

TPl«AkPH(3il)*TT<l) 
temp  f  0. 

Vl*ALPH(l,l)»OSpHl 

TlXiVX 

Vl?aAkPH(2,l)*DSDHl 
06  6o  Jal|4 

FAClaOELH( J)/(3,*6o76.1155) 

rAC2>2,«rACl 

IMAX*  t66{J) 

06  61  111, IMAX 
KaK*2 

HBH*DEkN(J) 

HlaH 

call  DDM(H) 

0RDH2f DRDh3 
DS0H2iDSDh3 
AN2aAN 
HaH*DELH(w) 

call  UDH(H) 

Y2*ALPH(1,K*1)*DS0H2 
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Y2ltY2 

Y22«ii PH(2iK*l)«0SnH? 

YJ«,4LPH(1,K,?)*DSDW3 

Y:2aAtPH(2iK*2)«DSDHS 

AN3*AN 

1»  (EL6v  ,LT.  1,  ,AND.  H  at.  201.)  5,  6 
9  CC«CS*CS*a,/RAD>HErO*QRAD/Rpl) 

CClal. /(CC*6076,li5b) 
pReo»<i 'Cc^Hi 

reiLewiNG  is  AppRsxtPATieN  required  near  tmetaio  anp  hiq  roR  range 

CALCUkATieN,  RANGE  IS  CALCULATED  THUS  F8R  H  •  lOO  AND  H«  200  WhEN 
elevation  angle  is  less  than  X  DECREE, 
RNQaCCl*PRSD/(S0RTr(PH8O«SS)*SN) 

DSlaRNC 

APPRexiHATE  ATTENUATIUN  IS  RANGE  (ThCoWAV)  TIMES  AVERAGE  VALUE  8F 

gamma  In  the  range  Interval, 

ATTENXaRNO*(ALPH(l,X)  *  ALPH(X,2)) 

ATTEN2aRNG«( ALRH(2,X)  ♦  ALRH(2,2I) 

tadar  Range  is  geehetric  range  times  average  refractive  index, 

RNGaRNGalRPl  •  X,  *  REFO aEXP ( •QR AOaHl ) )  a 

lSIGa2 

68  T8  7 

6  DSaRGKDRDHS) 

RNOaRNG*DS 

OSta  0$/( (ANt*AN2)*,5> 

AITENi«ATTENl*ATTl(Y3) 

YX«Y12 

Y2aY22 

ATTENJaATTENE  •  ATTX(Y52) 

7  RG<K*i)aRNG 

attn<x»k*x>«  atteni 

ATTN(2*K^;)PATTEN2 
ATTnHiM*X>*ATTENx  a  ATTEN2 
G9  Te  <10|li)  ISIG 
IX  PR8D*2>*CC*H 

RNOaCCXaRRSO/lSORTFIPRRDaSSIaSN) 

DS29RN0«0SX 

ATT6N;aRNG*(ALPM(l,i)  a  ALPH(X|3I) 

ATTeN;<aRNGa(4LRH«2.i)  •  ALRH(2,3I) 

RNOpRNCaIRPi  ♦  I,  a  REPoaEXP ( aORADaH  ))  a  ,9 

ISIG*1 
08  re  X2 

XC  0SaR02(DR0H3) 

RNQaRNOaDS 

DS2ffD8/( <AN2aAN9la.9) 

YXaYlX 

Y2aY2l 

ATTENl*  ATTEnX  ♦  ATT2(Y3) 

YXaYX2 

Y2aY22 

A1TEN2a  ATTEN2  •  ATr2(Y32) 

X2  RG(Ka2)aRNC 

ATrN(l,Ka2)9ATTENl 
ATTN(2iKa2)aATTEN2 
ATTN(9,Ka2)aATTENx  *  ATTEN2 
ALeSS»XO.**I'ATTN(3,Ka2>a,09» 


o  o  o 
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TP3aAt>^Ht3tK*2)*TT(«(*2)«tt.«SS 

OSl  •  (RG(K*1)*H0(K)  )/((ANi*AN2)«.)) 

OS?  ■  (RG(K*2)»R6tK*l))/(UN2*ANS)*,5) 

If  (K.EO.X)  70. T1 

•PPReXlMATlON  EMPLOYEU  in  place  SF  first  INTFGRATIRN  step, 

IflED  Te  Olvf  VALID  RESULTS  IN  MIGH.ATTENUAT IRN  CASES,  ANAL^TtC 
APPRSAIMATISN  STARTS  AT  STATEMENT  70, 

70  CEXfO,S*C«NST*(ALPM<S,l)«ALPH(3,2)) 
CEYaO,S«C|NST*(ALPH(5.»)*AlPHI3,3)) 

ALSSlfExPF (•CEXPDSl) 

ALSS2aExPF(*CEYaOSl) 

ALSS3|EXPF (•CEY«(DSX*US2)) 

DTEMP*<0.»/C«NST)*( (TT<x)*TT(2>)«(l,.ALeSl>  *  < TT { ? ) .TT ( S > ) • ( ALBS? 
1  *Ai,esS)) 
oe  TO  72 

71  ALSSSllO.**(«ATTN(3,K*l)«.OS) 
rP2iAtPH(3,K*l)«TT(K*l)*AL8SS 

CALL  JNTaRT(0Sl,DS2,TPl,TP2,TP3,0TEMP) 

72  TfeMP  I  TEmPaDTEHP 
DROHliORDHS 

TPl  t  TP3 
Yl«Yll»V3 

Y12«vA2 
ANl  ■  AN3 
ei  CONTINUE 
*0  CONTINUE 

ATMPaTEMPpCONST 

END 
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SUBRfUTINfe  ODH(N) 

C«nh«N  /RRQ/RErOiRiOtQ*AD>U 
crnmrn/drs/  dsonsiDrdhSi an 
EXiREro*EXP(>6RAD*N> 

AN«l,*Ex 

V*EXt(2.*EX) 

Wi«H/RAD 

W«W1*(2.*N1) 

0SDH3|AN*Ui*N1}/S0RTK  (U*V4H«V«H) 

DNDH3|AN*OS6h3 

END 
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SLHRCLTlNfc 

Ce^HPN  /PTe/  PP<75),  TT(7b),  RP ( 7i > , AlPH « J . 79 ) 

^6MMP^ /H20/RM€rAC 

Cl  PENSION  R  TRP(23)  iPTrk (33) ,nbLR(A  ) . IMAx(4) 
CATA(CELH«100..1''OC,i200(J.,iOOf',).(|f‘AX«?1.28.?Oi4> 
rATA<Ai.2»,l)i(Al.3i.27),(hlt260f’fl,).(W2»d?pno,)i(“RR=,5> 
rATA(CSTH2G,  A, 696  E e I ) , ( P fi-2Pt 22 . 236 > , ( UEl 2£Pe*l7 . 99  ) 

CA  T  A<  CFNSTt  3,7  i4P».(CSAST2io,206''£),(CV?»,  750044)  I  (S^-XI*.  20046) 

TA TA(f T RPa* 6, 264 "*50, 4466, 59. 59l ",60. ABAfljOi.i'iOAiAl. 0002,62.41121 
162 ,9960,63.6605,64,1272,64.4779,^5,3240,66, 7626.66,2978,66,031?, 
267,3627,67,8923,60.4205 ,60.9478,69 ,4741, 70, 0000. 70. 5249, 71.0497) 
TATaO  TuMi 11 6,76 06, 62. 4063, 6  0, 30  61,69, 1642,58. 32 39, 57, 6195,56,90 32 

I, 56,3634,65 ,7439,65, 221  4, 54. 6720  ,6 4, 1294. 63. ‘>960.53,0695, 59.5450. 
252,OP69,6l,5;4l.5O.9949.6O.483".49,973O,49.464«.4B.9502.40,453O) 

CaTa  (PP= 

II. ni326R‘3,i,C0969E*2,1.00596£»3,l.0Ot3lP*3,9,9a4,8i5E‘2,9,95076R«P, 
l9,9i473t‘2.9.67880e*2,9,64999E«2.5.6n72B£*2.9.77167E»2,5.7?6i7F»9, 
l9,7op7  7e*2,9,66540L*2.'5.6;<O2  9fc«2,9,5952iP«2,9.'..6O?3fc*2,9,59536R*2, 
i9.49ri5Re*i,9,4559iE.2,9,42i35e*2.9,oPl3o£*?.P.75l29E«2,0.43in9C«2, 
l6,:2047fc*2,7.ei9?ir*2,7.527i"t.2.7.2439iP*9,6.96943E.2.6. 7"347E.9, 
l*,A46aie*2,<.,l9629E*2,5,95459£*2.5.7?065R*2.5,4942?E*2.6,275l3P*?. 
l',O63l9fc.2,4,e5092F42,4.66og?fc«p,4,46(s46F*2.4.283i4F*2.«.iO449F»?, 
I3.93i76fc«2,3  ,  7649  7E«2  .  7  ,  6  0398fc*2.3.44662E.2,?.p9674E*2,3,l542  0F.2, 
l3,Pl4e4fc*2,2,75i3uF»2,;.50641t.2,2.27969F,2,?.07348E»2,l.e3230E.2. 
U,7io<3£*2,i,55420E^2,l ,4l94iE,2.1.2PJ52F,?,i  iaft4iE*2,l.06on0P*?. 

l9,633:9t.i,H,  75472r*l,  7.96650£*1.7,23n7F.i,6,S79i?F*i,  5.9  73;)3E»i, 
l5.429Cl£«l,  4,93  4  4  7F»l,4,4HSo6E*l,3,6s3282F*l,2.7a'(0  7E.l,  2,1949ifc*i, 
11, 7376 5fc*l,l, 3 »270F»i,i, 105338 *1) 

CATa  (7Tt 

l2,38l6nfc«2i2,e7962F»7,2 , 3 77e4E *2 ,i , e 76666 *2 , ? . ^ 73688 ♦£ . 2 .671690 *9 , 
l2.8697iE»2,2,e67  73E*2 .2 .P65  75e*2.<.8/37/E»2.?.H4j  79E*2.2.e598iP*2. 
19.857£3e«?,2.t55"9(*2.2,e6367E»2.2.65i8aF.j,?.h49OoE*2,2.04792E*9, 
19,R4594e.2,2,64396E*2.2,84«9«E*2,2.892l7F*2,7.fio?>>7E*2,2,7«256F*2, 
l2,762/6E.2,2.74296E*2,2,72tiA,E.2,2.7C337F*2,?.68357E*2,2,66378F*9, 
12. 64 3 99e *2,2. 6249 OE *2, 2. 60442E *2, 2. 5 8463F *9,2. 664 fl6fc*i, 2, 5<607F*?, 
l2.625  29E.^,2,6r|6‘'lF»2,;,4C5  74E.2.2.4  659;F.2,2,44620E*2,2.<2643F*9, 
l2,AO090E*2i2,3P0''7P*<«F.397i3£*2,2.34737R*3,?.3i761E‘4,2,3O785E*2, 
12,2e509E*2,2,24669E»2,<,20909E*2.2,io96'jc«?,2.5666Lc*2,2,16o60O«2, 
12. 16669R*2*2,166‘0E ‘2.2.104400 *2, 2,1* 06 0F*?,216660E*2, 2. 1*660F*2. 
12,366tnL‘2.?,1666uF*2,?,1686"c*2.2.a*66oF*?,2.1666"E*2.<,166600«P, 
12,16660E»2,2,16660E*2,2,l6*6"E«2,2.l4.66CF*?,2.366*OE‘2.2.i9069R»2, 
lP.2360PE*2i2,i8l''4P‘<.2,32664E*2) 
rATA(R«*  7,5CC00f-  0  .  7,4358*E  r. ,  7, 371790  Q.  7.30779E  ", 


1 

7.2438/6 

0  . 

7.10o:4E 

0, 

7.1343CO 

0. 

7,(152070 

0. 

0.9P935O 

0  , 

1 

6 ,9?67<jF 

0, 

6  .Pe?47E 

0. 

6,799330 

0. 

6, 73*336 

0, 

6.6734/0 

1 

6.61C77E 

0. 

6.548224 

C, 

4 . 4PSP5E 

0. 

6,42304E 

0. 

5.3*1*^ 

0  . 

1 

6, 299/66 

0. 

6.238144 

0  , 

5.633‘-fc£ 

0. 

5,  O544OF 

t). 

4.5.r6i5F 

0  . 

1 

3.989366 

0, 

3,507C2E 

0, 

3.0467.:F 

0. 

2,620394 

2^2* 6 20 r 

0. 

1 

1 . 90459F 

0. 

1.P1971E 

0, 

1,3'’C°£F 

0, 

1.1812"6 

0  r^l  .  0 1 0  9 1 F 

0  . 

1 

e , 646354 

•  5  . 

7,304476 

-1 . 

6, 291704 

-1. 

8,3431  3  6  ^-1, 

4.57l"6r 

-1 . 

1 

3,83698F 

*  1  . 

:.P9824L 

-1  . 

2 , 7A74tF 

-1. 

5.3i43(Hr 

*1  . 

1.93  73  7E 

-  1  . 

1 

1,560236 

*  1  . 

1  .734.J7E 

-1  . 

9 , 465r2f 

-2, 

7, 08o406 

6,233'5E 

•  2. 

1 

2  ,  eP4,i9F 

•  2. 

3 , 7u72CE 

-2. 

1  1o39''1F 

-2. 

6.453  99e 

.3, 

4, 012*40 

-3, 

1 

2,3P249F 

•  3  . 

1  ,466  706 

•7, 

1  .OSpPic 

-3, 

9,027o7E 

.4. 

8,406426 

-4, 

1 

7,915  78  0- 

-A. 

7  ,  3  <(4oC4 

'4  , 

6,376*86 

•4, 

8,799026 

-4. 

5.5)1070 

•4  , 
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1  5,4«J62E  .4,  5,?>UP5«£  --i,  •j.640?<fc  S,8797?E  -4.  6.2M32P  -4, 

1  7,067266  *4.  e.iJft4?6  '4.  V.C<»1‘'56  -4i  7,j2874E  -4,  ‘>,6'?1?26.  -4, 
I  4,4')<346  -4) 
r;AT*mMof  ic.i ,  > 

rGH2tf  ‘<H2«l  ,c.3 

TRATIGif  r.r1?/ffiH2'' 
fSiJM2*(rGH2*FNk'2''')*»? 
rcir2»trGH2.r4K2?j«-2 
1*0 

^•*•100. 

:e  loc  j*i.4 
I  P  *  I  f  A  X  (  J  ) 
fg  100  4*1. !>' 

I  *I  ♦! 

t-G  *30(1)  .  bh'>  AC 
pp-iRc  •! r  <  I ) /7ap , 7i» 

PPapp  (  ]  I  «  PP»,/CVP 
TbOtTT  (  n«*2 
p»p*rfcLH ( j 1 
If  (I-.L6.H1)  10.11 

10  Atl* . 6 4 

re  Tt  15 

11  If  (H  .3T,  P.?)  i7.1A 
1?  ALlsl . 157 

GG  Tt  15 

13  4ll*,s4  ♦  ,  7l  7  •  t  •'-P  i  ) /'4O0"  . 

15  P4LlsC.lVbT?»*Ll«'’''  (  I  )/TT  (  1  ) 

►•5Al*PAtl«^Al.i 
fo*>"*ui/<^rH’i«M^Ai ) 

«cf'*a , 

re  50  H,i,;j 

AM  »  AN  »1  , 

LNP2iA?v«(2,*AN»,1,  )/ANi 
l  (.►‘2  S  ANl  •  (  7  ,  •  A»  -1  .  )  /  AN 

C\n252,«(AA»i.\*A'*l,)»(2.»AN*l.)/(AN*ANl) 

fNPjPALl/ttfTr.Pv'-l-I.CM? 

rNPiP-Al,1/((pT«‘‘(‘'l.frM7  )»«2*P':Al)*pAll/f(f'-'“(P)»Rr.P2)»»2»nSAl> 

Tl.Rv:(F-,lP«i.  \iC2*f''''*C''r*7*f&«l,‘‘07>»t»pr  (•2.r,6tj44»A''*AM  /  TT  (  [  1  ) 
5LMfSC''*TfcK.'i 
50  '■GNTlN'Jt 

alPpG2  s  ,'tT  (  p..3)*rGh2^*^J“ 

ALPp  (III  )  *AL  =  pe2 
T40tJOn./T1  (  I  ) 

''fcLMCfcL/fe‘=J*tPP»»T'-!B«  PP(I)«CVP  •S4y]»Tw:-'»»,6;M»l,t-0 

OELf  ?»nkLf 

rPP2:(WATi;*:£Lf*<l./(f5>Cf<R*'’felf7)«i,/(rrjr2»DLLf2>) 
ALPH22erb7P2"*fO“^2»PP«*T‘’e**3.6*ixPM2.l44»(l.*TR.0>)»ft’W7 
*Lf  P6S*l  ,0/ i^^-?  •  f  '.P2?  •  P  -  •  PP  •  TT  (  I  )  .  .  i -2 .5  ) 

ALPw(2.|)  *  AlPh?A  *  AL^Pt'^ 

r  ai.ph(i,[)  15  AOS,  r/Rff.  f'^o  c^Y^tN,  Ai5p(?,n  katpr  vapbr. 

C  ALPh(J.I)  IS  total  ARSrKPTTJ.v  rcPFFirjfc,jY 
AlPh(3.1)»  ALPp(1<I)  ♦  *LPk(?.I) 

100  GGMIN'Je 
f  MI 
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SOBReu^INfe  |»iTGRT  ( HI .  U2 .  Y 1 1  V2  .  T3 ,  ARE  A  ) 

h22|H2«H2 
HHIHI aKZ 
HPHtHi*H2 

*f4r«(Yl«H2  *  Y3«m1»V2«HPH>/<HM«HPH) 

AREA  e  (Af AC/3.  )«tH22«H2  *  Hi2«Ml)  •  ( y3- V2- AF AC«H?2 ) * ( H?2-H12 ) / 
1  (2|*H2)  •  V2«hPH 

ENO 


o  o  o  o  o  o 
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SUBROU^INfe  IT6RaT(Rmx,4TT) 

GIVEN  Radar  maXJHuh  RANQE  |N  NflNABSQRetNG  SPACE,  RHXi  THE 

subroutine  finds  TNE  atmospheric  attenuation  ATT,  FROM  A  TABLE 
supplied  NT  anotmer  Subroutine  via  the  corro^  statement, 

IT  then  CBRRECTS  the  range  RHX  by  a  factor  based  bn  ATT,  THIS 
PROCEUURE  IS  iterated  UNTIl  SUCCESSIVE  ATT  VALUES  DIFFER  BY  LESS 
than  9,X  decibel, 

CBMMON/RGA/  RG(75),  ATTN(3,TS) 

ATTL«8, 

1  CB  XO  I>2,r9 

IF  {RS<M  •  PM*)  xo.o.ll 
iX  J«l»l 

ATT|(ATTN|3, I )-*TTN(S, J) )«(fiMX«R6{ J) )/(RCM )-RG( J) >  *  ATTNf3,J) 
GB  TO  12 
9  ATT|ATTN(3, I  ) 

GB  TO  12 

10  continue 

ATT,ATTN(3,75) 

X?  OIFF  I  ATTL  '  ATT 

RMX  •  RMX  •  10.«*{81Ff *,025 J 
If  (AHSF(DIFF)  ,i.T.  0.1)  RETURN 
IF  (RHX  ,6E,  HG|79))  RETURN 
ATTL  ■  m 
GB  Te  1 
END 
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1 7. 1 5( 1 


1 

2 

3 

20 


subroutine  POSN(POTiPF*iNPULSiKASE*SOB) 

niMENSlON  OBO<5) »ScOPE(S) 

COMMON  /POS/  PFl*ST tNLBST #kSUAST 
n4T*  %fL»ST-0.) .  (NLAST«0)  .  (KSLRST—l) 

1  f pnrAC«4«Bf 20 • fSO* t 13« • 13* )  ^  9 

’oaTA<OBMIN—30.)  t<OBMAXP50.»  •<OBl«0.».(0e2«fc.> 

IF  (POT.EC.POLASn  1»20 
IF  (PFA.EC.PFUAST)  2.20 
tr  (^PULS,EQ•^l.*ST)  3»20 
IF  iMSE.EO.ItSLAST)  RETURN 
POLASTpPOT 
PFL*STpPF* 

NLASTpNP'JLS 
KSLASTpKASE 
KpKASE  ♦  T 

CBl-OR?<K?-StOPE(K)*AL06lO(PULS».<POT-.5)-PoFAC(K)MPFA-8.)..A. 

?S”iCe«(OBMlN.OBMAX.Oei.OB2.4,l5.N01.SOB.POi,POT,PO» 

END 


1. 
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function  PO(SNOB) 

COM«CN/POS/FA»NiK*SE 

NP«N 

r4N«F* 

KASbMASC 

CALL  MAPSnR  (SNOBfNPtFANyKASfPOl) 

P0«P01 

CNO 


11  w  uu  o 


-4 
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SURBOUTINC  INVeBS(XMIN»XMAX*XL0«X»>l«NSta*Lt<4«N0Tt)ltFltFT«F) 
this  is  a  MOOirtCATIOK  (ABRIL  1«70>  OF  A  BRevlOuS  SUBROUTINE  NAMEO 
invert.  This  version  has  AOOITIONaL  RaRAMCTERS.  L*  V.  blare*  NRL. 
This  subroutine  finds  value  of  X  Thai  results  in  F(X)  ■  ft*  by 
ITfRaTICN  SaSEO  On  linear  INTERBOlaTiON/ERTRaROLATjON  from  BBevIOUS 
TVO  trials,  function  f  must  be  nonoionic. 

TEST  a  10.«*«-NS10) 

Fo  ■  FT 

IF  (FT  .EC.  0.)  Fo  ■  I. 

NO  I  a  I 

dELTa  ■  XhI  •  XLO 

Xl»XLO 

*?»XHI 

FlaF(Kl» 

F?»t  : *2) 

SLOPEa(F?-Fl)/(K<?-Xl» 

IF  (Slope  ,eu.  o.i  10.21 

10  ►  '‘AXaF  (XHAXl 
i  -IN«.'  (XMIN) 

^uOPE  ■  (F-IAX-FhIN)  /  (XHAX-XNlN) 

21  t  -*(;r  J-FT1»SL0PE  .Ll.  0.)  22i23 
2? 

ri  iF? 

X2«X2*0ELTA 

IF  (X?  .QT.  XBAXt  X2«XMAX 
F2aF(K2) 

00  TC  >l 

23  !F((FT.Fi1«SL0Pe  .UT.  0.)  2At25 
2*  x2>XI 

XlaXUOCLTA 

IF  <Xi  .LT.  XHIN)  Xl«XHIN 
F2aFl 
FlaF'xU 
oO  TO 
2S  *A«X1 
XRax2 
FAafl 
FRaF2 

IF  (ABSF(F2-FT)  .LT.  A8SF(FI«FTI)  7»B 
7  F22aF? 

F2aFl 

F1-F22 

x22»x? 

X?aXl 
XlaX?2 
<I0  TC  6 

1  FlaF(X) 
xlax 

TESTl  a  *bS( (F1-FT)/F0) 

IF  (TeSTI  ,LE.  TEST)  2.6 

2  BETuRn 

6  IF  (ROI  .CE*  LIm)  12*13 
12  PBinT  AO 

PBINT  41 .  LIM 

PBINT  42*  XMIN.  XMAXf  XLO*  XHlt  NSIG*  LIM*  NOI*  X*  FI,  ET 
HETURn 
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15  IP  (PI  .EC.  P2)  14*14 
14  IP  (Pi  .EC,  PB)  17.18 
IT  *e«Kl 

VpaiR 
80  TC  14 
18  X*all 
X2*S4 

18  Ka(X4*XBI«.4 

«0  TC  1 

16  xa(Xl-X2)*(PT-P2)/(Pl-F2)  *  X2 

IP  (X  .LT.  XA)  XaXA 

IP  (X  .07,  Xa)  XaXB 

NO I  a  NO I  •  1 

F2arl 

Xpaxl 

00  TO  I 

40  pokhat  (//•  message  prom  subroutine  INVERS  —  •//» 

41  format  (a  PlNCTION  INVERSION  n6t  ACCOMRuIShEO  wITmIN  SPECIPIEO 

I  *13.  •  iterations.*/) 

42  format  (/a  INVERS  PARAMETERS  MERE  a.AlElO.S.Px) *3(I3*2X) •e10.3* 

1  2(2x.El0.3)  //) 


non  nf%n  nt%n  n  nr%  nr* 
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e 

e 

c 

c 

c 

c 

c 

c 

c 

c 

e 

c 

e 

e 

e 


suB)*ourt*te  ma*siki  (SMOB**itrA«iuu«Aiii 


inm)Ts  ««e  ••  Sum*  sibnal-to^isc  »o«cm  baTio  ia  dccIkis  >• 

Nt  WU«MB  OB  ButSCi  IBTCONATCO  •• 

Ba.  BaLSC  ALABm  MOBMILITV.  CAMfOtCO  At  AOtOCiitl  VALUC  or  •0«CB 
OB  rcN  (C.o.t  ba  ■  ••  mcam  io**«i«o.i  Balm  alam*  mobboilitv  •• 

KAtCt  S*C«L1M  BkUCTWATlOB  BOMit  AltM  RAM  ■  0  BOB  NOBBLUCTUATInO 
OilTBuT  Br  It  BBOBaOILITT  OB  BtUCtlOR 

BaMO  OM  BBBBBar  rBITTCR  aT^  aBBlICO  BMTtICt  LBOOBAteBT,  NAhZO 
SURBfvTlRC  BABCU".  MOOtBUb  At  IBk.  Bt  L.  BLaBC.  »•»!$  MBtlON 
OAtCO  ABBlt  l«ri 

aBl  VCBtien  OCBIMO  Br  At  ArLM  ALRta  BMBBIB  IUrBCUM  COKIBti. 

RBL  BOO  CRAMBIO  TNlt  TO  BALM  ALABb  BBOBaBIlITt  lAt  OTBIMO  ABOVC) 
tOM  OtwCB  CMAtMCf  ALtO. 


ftOuBLf  BBCeltlOR  C*M,*aBa,AABBB*B«B«B«*A*rAkfl.CltttCBtTt 
»OUBLC  BBtcftlOk  OOAB*  OCBAL*  tb^BO.  .  ArB.  tR 

COMBlTC  RABCUBaBRC^lBB  MftCtlOR  BBO^ARUimt 

COWVCBT  t/«  lA  M  to  MMBICAL  ty%  Aa'IA 

tAiB  •  lA.««<tROB*.|i 

to  CCRyfAT  tBlt  tWBB^tlRl  *0  1  €*  tA  At  t«A  RaACU**  BALM* 

ALaBB  RURRfBr  0>AR«l  tii|  AfA*  ftAttRBR*  ?$  BEaO  •’  «OnC  •  # 


••000200 

•0000)03 
000004 30 
oooootoo 


•0000000 

•oooieoo 


•00*4 J 

IB  MOoe  IS  1*  CORvfBt  $t  r«  rBa*  UBOMAt  «a  BalM^ALaBm 
•BOfABlLltt  BAta(A  ri«A«  BABC<^  BALM*BLABB  VB^A 

tB  (BOOCI  000*  BOO.  ••• 

too  BAN  BOLOOlOCOLOOt.tl/tkOOtl.^MO-ek^Bt^BAi »r 
00  TO  oot 
Boo  ban  •  BA 

test  INBUTS 

•ot  TB(N)  00.00.2 
f  IB(BAIOO.OO.) 

)  TB(RASei  00.4.4 
4  IB(RAte.4|  t.s.oo 

eSTlBATC  BIAS  LCVCL 

s  ^NBR  •  0. 

A  COBB  ■  FAK 
AN  ■  N 
YABR  •  0. 

IF  (aPRCV  .EO.  N  .ANO.  BABBCt  .CQ.  Ba>  00  TO  TTt 
IFtN-lt)  7.T.0 

7  VBPB  ■  ENA(i.*2.2*CNBB/EN**( (O.O/S.Oi^.nitBeNBR) ) 

00  TO  11 

a  tOPB  •  EN»U.*1.3*ENBA/CN««(.S*.011»BNBHn 


••••noo 

••••1200 

••••1)00 

••••1400 

ooooitoo 

•0001000 

••••ITOO 

•OOOIBOO 

•OOOIOOO 

•0002000 

00002100 

•0002200 

•0002)00 

000S24O0 

00002800 

00002600 

00002700 

00002000 

00002900 

0000)000 

00003100 

0000)200 

0000)300 

0000)400 

O0003SO0 

00003600 


r»  n  o  (%  c% 
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COMPUTE  bias  LEveu 

11  €NP«  •  10.*«eNPB 
OAMPR  ■  OOAM(YBPR|N>1| 

PVP  •  ,9»»<l«/tMPR» 

8UML  •  SUl^LOOtM.l) 
tPlOAMPP.PYB)  10<l2tl2 

10  H  •  ,01 

00  TO  lA 

12  H  ■  .,01 

lA  ¥0  •  YBPR 

EO  ■  OCVAL  (YO'N.USUML) 

16  Yl  ■  yO*m 

El  •>  ntVAL<Yl,N-l,SOML) 

STEP  a  OAMpR  •  Ha(C0*ei>/2, 

IF(SISNP(l,,SlEP-PYe)*SIONr (l.,M))  IB.POdft 
IP  YO  a  Yl 
FO  ■  FI 
6AHPR  •  STEP 
<30  TC  16 

20  !F(M)  22t2Af2A 

22  YR  a  VI  -  MafPYP-STEP) /IOAMPM-STEPI 
(30  TO  30 

2A  YP  a  YO  ♦  Ha(PYB-GAMPR;/(STEP-QAMP«» 

30  RlAS  a  VB 
777  YB  a  BIAS 
NP9EV  a  N 
FaPREV  a  FA 

SELECT  N.S  CASE 

X  a  SnP 
K  a  R4SE«1 

<30  TO  (100«200«3O0iAO0«S00)  t  K 
CASE  0 

100  SUM  a  0, 

P  a  eN*X 

tF(YB>P«EX)  150,102,102 

102  K8  a  .(CN,l.)/2.  ♦  SQRTF( ( <EN.l.)/2,)aa2aPaYB) 
KS  a  XMAXQFtKSiO) 

AS  a  1 ,.GAM(Te,KS*N>liTN) 

T5  a  EVAL(PtKS)aOS 
0  a  GS 
K  a  KS 
TERM  a  TS 
TU  a  TN 

110  TEMP  a  SUP.TEPM 

IFtSUM-TEPP)  112tll6ill6 
112  SUM  a  TEMP 

1F(K)  116,116fllA 

llA  term  a  TEOMaFLOATF (K)a(0-TL)/«P»6l 
G  a  6-TL 
K  a  K.l 

TL  a  TUaFLOATF (k*N1 /Y0 


00003700 

00003000 

00003900 

OOOOAOOO 

OOOOAIOO 

O0OOA20O 

0000A300 

OOOOAAOO 

ooooasoo 

OOOOAAOO 

O0O0A700 

OOOOAOOO 

OOOOA900 

OOOOSOOO 

00009100 

00009200 

00009300 

000O9A00 

00009900 

00009A00 

00009700 

00009800 

00009900 

00006000 

00006100 


00006200 

00006300 

O0OO6AO0 

00006900 

ooecatoo 

00006770 

00006BOC 

00006900 

00007000 

00007100 

00007200 

00007300 

00007AOO 

00007900 

00007600 

00007700 

00007000 

00007900 

00008000 

00008100 

00008200 

00008300 

OOOOAAOO 

00008900 

00008600 

00008700 

00008800 

OOOOP900 


u  o  o 
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«0  TC  110 

no  TL  •  TN«YB/rL04Tr  (KS*K) 

K  •  BS*! 

0  ■  61*TL 

TCBM  a  T$«P*0/ (GS«rLO*TF <Kt ) 

1?0  TfMP  a  SU*'«TEt*M 

ir (5UP-TCPP1  IZZ.lOOtlOO 
l^t  SUM  a  TEMP 

TL  •  TL«YP/FLOaTF(R*N) 

K  a  P*l 

TfPM  a  TeRM*P«(0*TLI/(G«rcOATF (K)» 
rt  a  C*Tl. 

<10  TC  IZO 

ISO  KS  •  •!.  ■  tN/8,  *  $0nTF(EN»«2/4.*P*ve) 
as  apMtxOF (KS«0) 

OS  a  0AM ( V8tKS*N*l (TNI 
iriOS)  |74«iTa«iSS 
ISS  TS  •  CVAU (P*K$>aas 
0  a  GS 
TEPM  a  TS 
K  a  MS 
Tl  a  Tn 

1«0  TEMP  a  SUM*TCFM 

IKSLm-temPJ  )G2(166<166 
U2  SUM  a  T£mP 

TF(Rt  1«««164t16A 

164  TERM  a  TERM*rU0ATF(lO4(0»Tt)/(P»G) 

0  •  G*TL 

TU  •  TL4FL04TF (K*M«1)/Y8 
K  a  R.| 

OO  TC  160 

1G6  TL  •  TNavR/FlCAtF (KS«M 
K  a  KS*l 
0  a  GS-Tl 

term  a  TS*P*G/(0S«Fl04TF(K) ) 

ITO  temp  a  SUM  .  TERM 

IF<Slm,TEMP)  172,174,174 
172  SUM  a  TEMP 

Tl  ■  TLAYP/FLOATF (M»NI 

term  a  TERM»P«(o-Tl)/(G«FLOATF (M»ll  I 

0  a  G»TL 

K  a  M*] 

no  TC  170 
174  SUM  a  i,-5UM 
lOO  PM  a  Slla 
no  TC  40 

CASE  1 

200  TF(N-I)  210,210,220 
210  PM  a  FXPF(-Ye/<1,*X) ) 
no  TC  40 

220  temp  a  1,  *  l./(EM*X) 

PM  ■  1.  -  6AH(Ye.N-2fOUM)  ♦  E*PF«(eM»l,)«LOOF(TFMP).YR/(l,,ENa 
1  •GAM(Y8/7eMP,M«2,0UM) 

00  TC  40 


00004000 

00004100 

00004200 

00004300 

00004400 

00004SOO 

00009600 

00004700 

00004AOO 

00004400 

OOOlOflOO 

00010100 

00010200 

00010300 

00010400 

OOOIOSOO 

00010600 

00010700 

00010600 

00010400 

00011000 

00011100 

00011200 

00011300 

00011400 

OOOIISOO 

00011600 

00011700 

00011600 

00011400 

00012000 

00012100 

C0012200 

00012300 

00012400 

O0O12S0O 

OOO12GO0 

00012700 

00012600 

00012900 

00013000 

00013100 

00013200 

00013300 
00013400 
00013500 
00013600 
00013700 
00013600 
00013400 
OOnuoOO 
OTOIAIOO 
('0014200 
))  ^0014300 
00014400 
00014S09 


ft  nn  ft  ft  ft  rtnn 
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CISC  7 

900  iP(N>i>  jiotSiOtazo 
910  •  i*prj-78/n.*»i) 

00  TO  00 

920  ON  ■  ).  -  0AM(Ye/(l.*X) .N^ltOUMt 
00  TO  00 

CASe  3 

Aoo  tP(N»2)  «I0«A2o,A30 

*10  Oh  m  (l.»?.*X*T8/tX»2*)«*2l*C*PF(-2.»Tt/t*,««l » 
00  TO  OO 

*20  ON  •  (irSvB/(i.*x) )«exPF(«r8/iu*in 
00  TO  OO 

*30  C  ■  2,/<2.«Cn«X) 

O  •  l.«C 

Tr(rB*o*Ex)  **Oi*sOi*9o 

440  OUM  a  Oi 

▼eon  a  I, 

J  a  N 

*42  TEMP  a  SUfaTCAM 

iriSUM.TCMOi  444i4*6f446 
444  OUM  a  TEMP 

team  •  TEPMarBaO/ECOA^E ( J) 


00014600 

00014700 

00014000 

00014000 

OOOlfOOO 

00019100 

00019200 

00019900 

00019*00 

00019900 

00019600 

00019700 

00019000 

00019000 

00016000 

00016100 

00016200 

00016900 

00016*00 

00016900 

00016600 

00016700 

00016000 

00016000 

00017000 

00017100 


J  a 

OO  TC  **2 

*46  "N  a  1.  a  0*M(y8«Na2»CUMJ  ♦  ^•VB*ET4(.<TB»4a2> 

1  «  0*£v*i,<TP*N*l)*(l.*C*TB*<Et4-2»»*C/D)*9UP 

GO  TC  00 

490  an  a  I,  a  OAM(r8«Na9'0UM)  *  T8aEVAU(Te*Na3>aC/0 

1  ♦  EXPE«.C*tB-(EN-2.»*l.O0E(O»»«U.«C*TB-tEAa2.>*C/0» 

2  aO*M(V8*O.Na9,OOM) 

00  TO  *0 


900 


501 

502 


CASE  4 

SUM  a  0, 

C  •  2./«2.*X> 
n  a  l.ac 
<5  a  C/O 
P  a  C«T8 

*9  a  O.afN* (T0»0l 1 /2.«SO«Tf ( (EN-1 . • C V8*0»  > ••2/4*» •TBaO> • lEN* I , I  1 
as  a  XH|NOP(KSfN) 

NS  a  xmaxoE(XSiOI 
K  a  KS 
J  a  N-nS 
XKS  a  KS 
K  a  XMINOF(NS.N) 
tE(re-EN*<l.*0) )  S90tS01«901 
OS  a  1,  «  GAM(Pt2*N-l«NS»TN) 
lE(GS)  526tS26fS02 

TS  a  CXPF (FXS«LOGnc> ♦ <EN-FNS>aCOCE »OI aSUMLOOIN) aSUMLOGlKS) 

)  -SUFloO(J)*UOOF(GS)) 

0  a  CS 


173« 
00017400 
00017900 
00017600 
00017700 
00017800 
00017900 
00018000 
00018100 
00018200 
00018300 
00018*90 
00018500 
00018*00 
00018700 
00018800 
00018900 
00019000 
00019100 
00019200 
00019300 
00019400 
00019990 
00019600 
00019700 
00019800 
00019*00 
00020000 
00020100 


I 
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term  ■  TS 
TL  ■  TN 

510  temp  ■  SUm-TEHM 

IFtSCM.TEMP)  S12tS16iSl« 

512 

irtK)  51«.S16«SI« 

5l*  TL  •  TL»P/fLO*Tr (2»N-R) 

TERM  •  TEPM*rLO*Tr(K)«(Q.TLt/(a*fLOATf  <N.K*H«Oi 
0  a  0*T|. 

K  ■  K.l 

AO  TC  510 

5l6  TF(KS.N)  5l8tS2«.526 

51P  TEBM  a  TS«0«FL0*TF(N-KS»«(aSatN)/(FC0*TPlKSal)anS> 

A  a  CS«Tn 

Tl  a  TNari.0*TF(2»N-l-KS)/P 
K  a  K5*l 

520  temp  a  SUM*TE«M 

TP(SCm.TEmP)  S22>S26iS26 
522  5UM  a  TEMP 

TF(K«M  S24iS20«520 

52*  term  a  TEOMa<lari.O*TF<K-»<>»(0aTi.»/(FCO*Tr<Kal»«G) 

0  a  6«TL 

TL  a  TL»FL0*TF(2aN-l-K)/P 
K  a  a«l 
AO  TC  520 
52*  P»«  •  SUM 
AO  TC  0« 

550  «S  •  0*M(P,2**-i-KStT*» 
tP(8S)  S7*»57*t552 

552  TS  a  F *PF (f uSaLOUE <CI ♦ (CN-F*S) aUOOf lOI aSuatOO IN) .SUMLOb t«S) 
1  •$UM|.00(JI*k.00P(0S)) 

0  a  05 
TEPM  a  TS 
TL  ■  TN 

$1)0  TEMP  a  SUaaTEHM 

IF  (SL'M.TEMPt  S*2*56*ta66 
5*2  SUM  a  TEMP 

IF(«)  S*e«S*6«56* 

"  »  tl  a  TL*P/FL0*TF (2«NaK) 

TEBM  a  TeOMaFtO»TF(K»«»0-TLI/«Q«PLO*TF(N-Fal)aG) 
r,  a  C-TU 
K  a  K»1 
AO  TC  560 

56*  TF(KS-N1  a6aiST6<S76 

566  TEBM  a  TSaQaFLO*TF<N-KS)a(OS*TN»/(FtO*TF|KS*l>aoS» 

A  a  6S»TN 

Tl  a  TNaFLO*TF (2*N-1-*S) /P 
K  a  KS»l 

570  TEMP  a  SUMaTEMM 

IF(SUm.TEmP)  5T2i576f576 
5T2  SUM  a  TEMP 

1F(K-M  ST*«576,576 

57*  TEBM  a  TEBM»Q*FU0*TF(*-S)*(0«7Ll/(FLO*TFtKai)a0) 

A  a  OaTL 

Tl  a  TLaFLOATF (2#N-l-5)/P 
K  a  Mai 


00020200 

00020300 

00020*00 

00020500 

00020600 

00020700 

00020900 

00020000 

00021000 

OOOIllOO 

00021200 

00021300 

00021*00 

00021500 

00021600 

00021700 

00021600 

00021900 

00022000 

00022100 

00022200 

00022300 

00022*00 

00022500 

00022600 

00022700 

00022900 

00022000 

00023000 

00023100 

00023200 

00023300 

00023*00 

00023500 

00023*00 

00023700 

00023900 

00023900 

0002*000 

0002*100 

0002*200 

0002*300 

0002**00 

0002*500 

0002*600 

0002*700 

0002*600 

0002*900 

00025000 

00025100 

00025200 

00025300 

00025*00 

00025500 

00025600 

00025700 
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(10  TO  870 
876  PN  ■  1,-SUH 
(JO  TO  80 

8€T  PB0BAeiH7< 

40  tPiPM  8l»»**82 

41  PN  P  0. 

OO  TO  94 

42  1P(P6«1»>  44t94«93 

43  PN  •  It 

04  return 


C 

c 

c 


error  HCSS*0E  POR  8*0  INPUTS 


04  MRITE  (61t9l  N,F4  tSNRfKASE 

4  E0RM4T  <iNo  /Son  unreasonable  call  Sequence 

1  7hS  OIVEN  //4«  N  •  IBtSHtSHEA 

f  E16.S|5*i6«R*SE  •  I8t 

PN  ■  0. 

PUS  a  Ot 

return 

ENO 


00028BOO 

00028900 

0002*000 

00026100 

00026200 

00026300 

00026400 

00026800 

00026600 

00026700 

00026SOO 

00026900 

00027000 

00027100 

00027200 

00027300 

TO  MARCUM.  ZERO  RESULT00027400 
a  rl6.R«8X*SHSNR  a  00027800 
00027600 
00027700 
00027A00 
00027900 
00020000 
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rUNCltOM  OOAMIB.N) 

00028100 

nOUBLE  PRECISION  SUM.  TeRM.TtMR.FJ.OOtM.  nt'^klt  B.  SUML.SUMLOO 

00028200 

Integral  ■  i>(sum.  jbo  to  n«  op  cpPR(J*LOOF(Bt>B.LoaF(NPAc 

((<28300 

SUM  ■  0. 

< ( (28400 

K  ■  B 

( ( (28800 

IF(K«N)  100i200.200 

( ( (28000 

100 

J  ■  N«1 

(((20TOO 

SUML  ■  SUmlOO(O) 

( ( (28800 

TERM  •  QCwAL(B«J«SUML> 

(((28900 

10 

TEMP  •  SUM ♦term 

(((29000 

IFISUM'TEMP)  18. 20.20 

(((29100 

18 

SUM  •  TEMP 

( ( (29200 

J  ■  J»1 

( ( (29300 

FJ  •  J 

( ( (29400 

TERM  •  Te«M*B/FJ 

( ( (29500 

00  TC  10 

( ( (29000 

20 

OOAM  ■  SUM 

< ( (29700 

return 

( ( (29800 

200 

J  •  N 

(((29900 

SUML  •  SUMLOO(J) 

( ( (30000 

TERM  ■  0EVAL<B. J.SUML) 

( ( (30100 

30 

Temp  ■  SUm.TCRm 

(((30200 

tF(SUM. TEMPI  35.*0.*0 

( ( (30300 

35 

SUM  ■  TEMP 

( ( (3n400 

tF(Jol)  40.30«36 

( ( (30800 

3* 

FJ  •  J 

( ( (30000 

TERM  •  TEfiM«Fj/fi 

( ( (30700 

J  •  j.l 

( ( (30800 

00  TC  30 

( ( (30900 

*0 

OOAM  •  1,.SUM 

( ( (31000 

return 

(((31100 

ENO 

( ( (31200 
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fUNCTlON  OeVAL  (Y,M|SldML) 

OOUSIC  PReC!S10N  XPON«CN«oeVALt  YtSUML 
xPON  ■  •y 
rriM)  20«20fl0 
10  fN  •  ^ 

xPON  ■  XPCN«EN*OI.OO<Y)«SUNL 
20  OfVAL  •  OEXP(*PON) 

PC turn 
CNO 


I  01300 
(  01*00 
I  OlSOO 
(  OlOOO 
(  OlTOO 
(  01800 
I  01900 
(  02000 
(  02100 
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function  C-*M(eiNtTN, 

single  precision  version  op  CG»h 

SUM  •  n, 

X  ■  0 

TF(K-N)  100*200t200 
loo  I  •  N • 1 

TERM  ■  EVALOtJ) 

TN  ■  TERM«FuOPTF ( J) /b 

10  temp  •  Slm«I£Rm 

TF(Slm«TEmP1  1S»20>20 
15  SUM  a  TEMP 
J  a  J»1 
FJ  a  J 

term  ■  t£fiH«0/Fj 

GO  TC  10 
?n  gam  a  SUM 

return 

2flO  t  a  \ 

TERM  a  EV4L10*J) 

TN  a  term 

in  tfmP  a  SUM^TEniM 

IF  (SUM-TEmP)  35tA0«*O 
OS  SUM  »  TEMP 

!F(J«ll  An«36«36 

3N  FJ  a  J 

TERM  a  TERMaPj/B 
J  «  .-1 
GO  TC  30 

*0  GAM  a  l,,S.iM 

return 

FNO 


( ( (32200 
( (02300 
(  (02400 
:  (02500 
(  (02400 
(  (  02700 
( (  O2B00 
(  (  O2Q00 
(  (  03000 
( ( (33100 
(  (0  3200 
( ( (33300 
(  ( (33400 
(  (03500 
( ( (33400 
( ( (33700 
( ( (33R00 
(  ( (33900 
(  ( (34000 
(  ( (34100 
(  (04200 
(  (04300 
(  (04400 
(  (  (34500 
(  (  (34600 
(  (  (34700 
(  (  (34P00 
(  (  (34900 
( ( (35000 
( ( (35100 
( ( (35200 
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(  (05300 

FUNCTION  eVALlT.N) 

(( (35*00 

kpon  •  •y 

( (05500 

1F(N)  20t?0»10 

(( (35600 

10 

xPo"  a  xPCN*EN»LOGF(Y)-SOHlOQIN) 

(( (357C0 
( ( (35800 

?n 

FVAL  a  EXPF(XPON) 

( ( (35900 

return 

( ( (36000 

52 
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FUNCTION  SOMLOO(N) 

OOUBLE  P«EC1S10N  *«  Bi  SUMUOO 
OIMENStON  A(lUOO) 
n«T<(OUMA  ■  0>)«(OUMB  ■  0.) 
NNAXaJOOO 

IF (OU»a»OumB)  20«10«20 
10  nUMt  •  I. 

OUMS  •  0. 

N(.*ST  •  I 
4(1)  ■  0. 

20  NN  •  xAqSF(N) 

IF(N^.n  3(j--:-40 
30  SUHUCO  ■  0. 
hetupn 

*0  IF (NN-NLASTI  $0>S0i60 
so  SUHLC6  ■  4(NN) 
aCTURN 

60  F  •  NL4ST«1 

tF(NK.NM*,)  70,70,80 
TO  no  72  !»K,NN 

T2  4tl)  a  .  OLOGlFtOATFd)) 

NLAST  a  NN 
GO  TC  50 

80  TF(N|.AST-KMAX)  a2,90>80 

82  no  8*  IaK,NN4x 

84  4(1)  a  4(I~l)  *  OLOG(FLOaTF(I)  ) 
NLAST  a  NkAX 
90  8  a  4(NNA«) 

F  a  NwaX*1 

no  92  IbK.NN 

92  8  a  r  ♦  0).OG(FuOATF(I)  ) 

suMtCG  •  e 

8ETUfi^ 

ENO 


08/12/T2 


<  ( (36100 

<  (  (36200 

'  (  (36400 

<  ( (36600 
(  (  (36700 
(  (  (36800 
(  (  (36900 
(  <  (37000 
{ ((37100 
(  ( (37200 
(  (  (37300 
‘  < (37400 
(  ( (37500 
( ( (37600 
( ( (37700 
(  (  (37800 
(  (  (37900 
(  (  (38000 
( ( (38100 
( ( (38200 
( ( (38300 
( ( (38400 
( ( (38500 
( ( (38600 
( { (38700 
( ( (38800 
( ( (38900 
( ( (39000 
( ( (39100 
( ( (39200 

( ( (39300 
( ( (39400 


NRL  REPORT  7448 


53 


REFERENCES 

1.  Blake,  L.V.,  “A  Guide  to  Basic  Pulse-Radar  Maximum-Range  Calculation,”  Part  1, 

NRL  Report  6930,  and  Part  2,  NRL  Report  7010,  Dec.  1969. 

2.  Bean,  B.R.,  and  Dutton,  E.J.,  “Radio  Meieorology,”  National  Bureau  of  Standards 
Monograph  92,  U.S.  Government  Printing  Office,  Washington,  D.C.,  1966. 

3.  Blake  L.V.,  "Machine  Plotting  of  Radio/Radai  Vertical-Plane  Coverage  Diagrams,” 

NRL  Report  7098,  June  1970. 

4.  Boothe,  R.R.,  “A  Digital  Computer  Program  for  Determining  the  Performance  of  an 
Acquisition  Radar  through  Application  of  Radar  Detection  Probability  Theory,” 

U.S.  Army  Missile  Command  Report  No.  RD-TR-64-2,  Dec.  1964. 

5.  White,  D.M.,  “Radar  Simulation  and  Analysis  by  Digital  Computer,”  Johns  Hopkins 
University  Applied  Physics  Laboratory  Report  TG-952,  Jan.  1968. 

6.  Fehlner,  L.F.,  “Marcum  and  Swerling’s  Data  on  Target  Detection  by  a  Pulsed  Radar,” 
Johns  Hopkins  University  Applied  Physics  Laboratory  Report  TG-451,  July  1962. 

7.  Kirkwood,  P.K.,  “Radar  Cumulative  Detection  Probabilities  for  Radial  and  Nonradial 
Target  Approaches,”  Rand  Corporation  Report  RM-4643-PR,  Sept.  1965. 

8.  Nolen,  J.C.,  and  others,  “Statistics  of  Radar  Detection,”  Bendix  Corporation,  Baltimore, 
Md.,  Feb.  1966. 

9.  Killinger,  D.K.,  “Computer  Aided  Radar  Design  (CARD),”  Naval  Avionics  Facility 
Interim  Report  TR  1461,  12  Aug.  1969. 

10.  Minzner,  R.A.,  Ripley,  W.S.,  and  Condron,  T.P.,  “U.S.  Extension  to  the  ICAO  Standard 
Atmosphere,”  U.S.  Dept,  of  Commerce  Weather  Bureau  and  USAF  ARDC  Cambridge 
Research  Center,  Geophysics  Research  Directorate,  U.S.  Government  Printing  Office, 
Washington,  D.C.,  1958. 

11.  Sissenwine,  N.,  Grantham,  D.D.,  and  Salmela,  H.A.,  “Humidity  up  to  the  Mesopause,” 
Air  Force  Cambridge  Research  Laboratories  Report  AFCRL-68-0550  (Air  Force 
Surveys  in  Geophysics,  No.  206),  Oct.  1968. 

12.  Blake,  L.V.,  “A  Modified  Simpson’s  Rule  and  Fortran  Subroutine  for  Cumulative 
Numerical  Intc“gration  of  a  Function  Defined  by  Data  Points,”  NRL  Memorandum 
Report  2231,  Apr.  1971. 

13.  Meeks,  M.L.,  and  Lilley,  A.E.,  ‘  The  Microwave  Spectrum  of  Oxygen  in  the  Earth’s 
Atmosphere,”  J.  Geophys.  Res.  68  (No.  6);1683-1703  (Mar.  15,  1963). 

14.  Reber,  E.E.,  Mitchell,  R.L.,  and  Carter,  C.J.,  “Attenuation  of  the  5-mm  Wavelength 
Band  in  a  Variable  Atmosphere,”  IEEE  Transactions  AP-18  (No.  4);472-479  (July  1970) 

15.  Blake,  L.V.,  “Radio  Ray  (Radar)  Range-Height-Angle  Charts,”  NRL  Report  6650, 

Jan.  22,  1968,  and  Microwave  J.  11  (No.  10);49-53  (Oct.  1968). 

16.  Blake,  L.V.,  “Ray  Height  <”orr'.p'.'tafi''n  for  ?  Continuous  Nonlinear  Atmospheric 
Refractive-Index  Profile,”  Radio  Science  3  (new  series)  (No.  l):85-92  (Jan.  1968). 


Appendix  A 

A  FORM  FOR  LISTING  DATA  CARD  INPUTS  TO  PROGRAM  RGCALC 


The  following  form  has  been  devised  for  listing  the  input  quantities  needed  on  data 
cards  for  Program  RGCALC.  The  form  provides  a  convenient  reference  for  punching  the 
data  cards,  since  it  indicates  the  card  columns  for  each  input  quantity.  The  “case”  options 
(1  through  7)  are  also  described  in  a  footnote  on  the  form. 


DATA  CARD  FORMATS  FOR  PROGRAM  RGCALC 
(Reference:  NRL  Report  7448) 

CATdt  1,  3,  5,  7  .  .  :  Name  of  radar  (Format  10A8).  Cult  1-80.  <Any  other  deacnplive 
alphanumeric  matenal  cam  be  entered  in  theie  columna.) 

Cardi  2,  4,  6.  8  ...  ;  Radar  daU  aa  below  (Note:  A  quantity  with  F  format  can  be  punched 
without  decimal  point  if  it  has  an  integer  value  and  is  ngh‘  adjusted.  Non-integer  numbers 
must  be  punched  with  decimal  point:  right  adjustment  tht  unneceitary.  Integer-format 
quantities  must  be  right-adjusted  ir  specified  column  field  ) 


DATA  ITEM 

”1' 

...L 

FORMAT 

SPEC. 

COLS.  1 

Trinjmitur  power.  kW  fp,) 

F6,0 

Fulse  length.  ;j»oc  (r) 

—  .-4- 
1 

.  ,  .1. 

1 

F6.0 

7-U  J _ 

13-16  T 

Transmit  anienra  gam.  dB  ((7() 

F40 

Receive  antenna  gain.  JB  (C^) 

F4.0  '  17  20  1 

Target  cross  section,  sq  m  (o) 

F6.0 

21-26  i 

Frequency,  MHt  if) 

.  ...L 

F6  0  1  27-32  j 

Antenna  ohmic  lo«.  dB  (Lq) 

r 

F4  0 

33-36  1 

Receiving  line  lose.  dB  {L,} 

F40 

37  40 

Transmit  line  Inu,  dB  (Lt) 

1 

F40  '  41  44  j 

Antenna-pattern  scan  Iota,  dH  Up) 

F4  0 

45-48  ] 

Miscellaneous  loss.  dB  (Lg) 

F40 

49-ir!^ 

Bandwidth/ihape  factor,  dB  (C^) 

1 

F4  0  53  56  ! 

Receiver  noise  factor,  dB  iSF'. 

F4  0 

57  60 

Number  of  pulses  (M) 

1 

h 

See  Report  6930.  Eqi.  (71)  and  (72) 

..  J 

15 

61-65  ' 

Cases  0-5:  Probability  of  detection 

C'ise«6-7:  Signal-to-noise  ratio,  dB 

F4.0  66-69 

Cases  0-5'  False  alarm  exp't  (~logio 

Pfa) 

Case  6:  Blanlt,  Case  7,  Signai/noise,  dB 

F4  0 

70-73  j 

Case  (0  lo  7)* 

1 

11  ■  74  ; 

Target  elevation  angle,  degrees 

...  1 

F4  0 

75-78  1 

Calactir-noise-level  code  (“1,  0.  ♦1) 

I 

(minimum,  average,  maximum) 

.  1 

12 

RADAR  NAME  RADAR  NAME 


■I - 


H 


_ I. 


V*  0  lo  4  ‘  rang*  it  calculated  for  corretpondmf  S«erl.n|  nurtuation  rmmr.  S  -  calculated  for  all  fivr  Swerlmg  caa««  (0.  1 ,  2.  3,  4). 
S  “  calculated  for  the  S/N  value  in  Cnlt  7  -  calculated  for  Iwo  S/N  values,  on*  m  66  69.  other  in  70  73. 
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This  is  a  study  of  the  penetration  and  craciriog  of  rocks  using  lasers  available  at  the  Avco  Everett  Research 
Laboratory.  The  objectives  ate:  1)  to  obtain  data  on  the  rate  of  rock  damage  fot  various  laset  conditions  and 
2)  to  present  an  analytical  program  to  predict  the  temperature  and  stress  in  rocks  for  pointed  and  annular  beams 
of  radiation. 

The  laser  power  outputs  used  thus  far  was  fiom  I  to  IT  kW,  CW,  10.6  microns  and  pulsed  last,  power  up  to  1000 
joules  in  20  micioseconds.  Data  was  taken  with  sharply  focused  as  well  as  dcfocussed  beams.  Data  is  presented  for 
continuous  irradiation  as  well  as  pulsed,  pointed  beams  and  annular  radiation  patterns.  Three  type.s  of  hard  rock 
were  tested  namely;  quattzite,  a  Rhode  Island  granite,  and  Dresser  basalt.  Hole  penetration  energy  per  cm  of 
penetration  was  found  to  be  independent  of  laser  intensity  over  6  orders  of  magnitude.  The  specific  energy  for  rock 
removal  was  found  to  depend  on  the  fifth  root  of  laser  intensity  ovet  the  same  range. 

A  computer  program  was  developed  using  the  finite  element  method.  It  takes  a  radiation  pattern,  and  delivers 
curves  fot  tbermoclastic  strain  and  stress  in  any  direction.  This  program  is  capable  of  working  with  high  intensity 
laset  radiation  because  it  includes  nonlinear  heat  conductivity,  a  temperature  cut  -  off  at  the  boiling  point, 
anisotropic  elastic  conditions  and  a  special  iterative  procedure  to  avoid  computation  instabilities  resulting  from 
too  rapid  heating.  This  program  predicts  that  efficient  cracking  of  rocks  will  be  achieved  by  directing  the  radiation 
in  a  large  annulus.  The  beam  is  moved  so  that  the  temperature  of  any  point  stays  below  the  melting  point.  However, 
thermoelastic  stress  builds  up  in  depth. 

These  predictions  have  been  verified  by  experiments.  However,  the  experiments  to  spall  away  basalt  and 
quartzite  consume  more  energy  than  predicted  from  the  thermoelastic  computer  program.  This  is  believed  to  be 
because  the  computer  predicts  the  condition  of  crack  generation  whereas  the  experiments  show  the  complete  separation 
of  rock  material  when  the  cracks  are  extended  to  complete  spall. 
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